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On  this  grant  we  made  major  progress  in  three  areas: 

1.  Adiq)tive  G>ntrol  of  Nonlinear  Systems 

2.  Approximate  Linearization  (by  state  feedback)  of  nonlinear  systems 

3.  Software  tools  for  CAD  of  nonlinear  control 

Adaptive  Ctmtrol  of  Nonlinear  S3^tenis:  In  this  work,  we  extended  our  previous 
work  on  direct  adaptive  control  of  Single  Input  Single  (Xitput  nonlinear  systems  to 
schemes  for  adaptive  identification,  indirect  active  control  and  also  adaptive  model 
matching  of  Multi  Input  Multi  Output  nonlinear  systems.  We  also  studied  adaptive  ver¬ 
sions  of  the  nonlinear  regulator. 

Approximate  Linearization  (by  state  feedback)  of  nonlinear  systems:  While 
the  full  set  of  conditions  for  input-output  linearization  of  a  nonlinear  system  by  state 
feedback  have  been  given  in  the  literature,  the  question  of  how  to  proceed  when  the 
conditions  follow  slightly  short  of  being  met  have  not  been  answei^.  For  example, 
input-output  linearization  hinges  on  a  certain  set  of  regularity  conditions  (existence  of 
relative  degree  in  the  SISO  case)  and  minimum  phase  conditions  being  met  by  the 
plant.  If  the  plant  is  not  regular  and  is  slightly  nonminimum  phase  the  techniques  of 
input-output  linearization  need  to  be  modified.  We  discussed  these  techniques  in  the 
context  of  flight  control  and  also  other  examples,  for  instance,  the  ball  and  beam  sys¬ 
tem.  This  in  turn  led  to  a  deeper  understanding  of  the  structure  of  the  zero  dynamics 
of  a  nonlinear  system  and  their  structure  under  perturbation. 

CAD  tools  for  nonlinear  controller  design:  We  have  deveiop«l  a  set  of  CAD 
tools  for  linearization  and  approximate  linearization  of  nonlinear  systems  using  spline 
software  which  operates  in  real  time  and  is  capable  of  accepting  nonlinear  system 
description  in  numeric,  tabular  or  functional  form.  A  usev  interface  is  being  written  and 
it  is  being  tried  out  on  several  examples. 
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In  the  years  of  this  grant,  of  the  students  supported  on  the  grant  I  have 
had  two  Ph.  D.,  dissertations  (S.  Behtash  and  J.  Hauser),  two  M.  S.  Plan  II 
reports  (R.  Kadiyala  and  A.  Teel)  completed  and  three  Ph.  D.  dissertations 
are  nearing  completion  (R.  Kadiyala,  A.  Pradeep  and  A.  Teel,  all  expect  to 
graduate  between  January  and  June  1992).  The  grwit  has  enabled  the  PI  to 
switch  my  focus  from  linear  adaptive  control,  which  was  the  work  supported 
by  the  previous  ARO  grant  to  me,  ARO  DAAG  85-K-0572,  to  several  areas 
of  nonlinear  and  adaptive  control  and  more  recently  to  the  development  of 
CAD  tools  for  nonlinear  control  systems  design,  a  project  which  continues 
with  the  newest  ARO  grant  DAAL  03-91-G-0171.  We  have  made  several 
trips  to  the  US  .Army  Ordinance  Research  Center  at  Picatinny  Arsennal,  NJ 
and  have  set  up  a  good  working  relationship  with  the  group  of  Dr.  Norman 
Coleman  and  a  design  project  in  fire  control  for  Apache  helicopters  which  is 
supported  at  Integrated  Systems  Inc.,  Santa  Clara,  California. 

The  work  done  on  this  grant  has  had  major  impact  on  the  field  in  two 
areas:  the  adaptive  linearization  of  nonlinear  systems  (publications  [2]  and  [9] 
of  the  list  below)  have  begun  a  new  field  of  research  which  has  advanced  the 
theory  of  nonlinear  control  and  has  important  implications  for  CAD  tools 
for  nonlinear  control  systems  design;  and  the  approximate  linearization  of 
nonlinear  systems  (publications  [1],  [3],  [10]  and  [11])  has  opened  new  lines  of 
investigation  in  developing  a  nonlinear  control  systems  design  methodology. 
Finally,  we  have  begun  developing  software  tools  for  CAD  of  nonlinear  control 
system.’?  (publications  [12, 13])  below.  A  brief  outline  of  the  areas  of  research 
fir  dint's  is  now  given: 
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Brief  Outline  of  Research  Findings 

1.1  Adaptive  Control  of  Nonlinear  Systems 

There  has  been  a  great  deal  of  recent  progress  in  the  development  of  basic 
theory  for  the  input  output  linearization  of  a  large  class  of  nonlinear  systems 
by  state  feedback.:  The  chief  drawback  of  these  techniques  has  been  that 
they  rely  on  the  e.xact  cancellation  of  nonlinear  terms.  When  the  nonlinear¬ 
ities  are  not  known  exactly  adaptive  control  may  be  used  to  asymptotically 
make  the  cancellation  exact.  In  past  work  supported  by  the  grant  we  had 
developed  a  direct  adaptive  control  algorithm  for  this  purpose.  In  work  in 
this  time  period  we  have  proposed  a  variety  of  other  schemes,  referred  to  as 
indirect  and  semi-indirect  adaptive  control.  In  the  former  scheme,  we  de¬ 
veloped  a  number  of  different  identification  techniques  for  nonlinear  systems 
and  coupled  them  with  the  input  output  linearizing  control  law  using  the 
certainty  equivalence  principle.  We  gave  conditions  for  the  convergence  of 
the  scheme  and  showed  that  it  had  several  important  advantages  over  the 
direct  scheme  including  no  need  for  over-parameterization.  The  conditions 
for  convergence  were,  however,  far  more  restrictive  than  those  for  the  direct 
scheme.  Consequently,  wc  also  proposed  a  semi-indirect  scheme  which  com¬ 
bined  several  of  the  attractive  features  of  the  direct  and  indirect  schemes  and 
gave  a  convergence  proof.  We  have  been  comparing  all  of  these  schemes  with 
nonadaptive  alternatives  such  as  sliding  i.node  control  on  several  examples 
such  as  induction  motors,  (publications  [6],  (12]) 

We  also  began  the  study  of  adaptive  control  of  MIMO  adaptive  nonlinear 
systems  :  in  this  research,  we  (joint  work  with  M.  Di  Benedetto  from  the 
Universita  di  Roma)  studied  two  schemes  for  the  adaptive  tracking  control 
of  MIMO  systems  with  parametric  uncertainty  in  their  dynamics.  The  first 
approach  is  an  adaptive  version  of  a  static  feedback  law  for  tracking  control 
based  on  some  results  on  asymptotic  model  matching  recently  proposed  by 
Di  Benedetto.  This  scheme  is  based  on  some  some  results  on  asymptotic 
model  matching  recently  proposed  by  Di  Benedetto.  This  scheme  is  based 
on  some  new  techniques  for  extending  the  so-called  zero  dynamics  algorithm 
of  Isidori  and  Byrnes  to  problems  of  stable  model  matching  followed  by  their 
specialization  to  tracking.  The  second  scheme  is  an  adaptive  version  of  a 
dynamic  precompensation  law  of  Descusse  and  Moog  for  linearization  using 
dynamic  state  feedback,  (publication  [9]  below). 
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The  schemes  are  in  the  spirit  of  our  earlier  work  on  adaptive  linearization 
of  nonlinear  systems,  done  on  this  grant  with  Isidori  and  Kokotovic,  (publi¬ 
cation  [1]  )  which  were  however  confined  to  the  SISO  case.  These  schemes 
represent  the  initiation  of  a  full  theory  of  Model  Reference  .Adaptive  Control 
of  MIMO  nonlinear  systems.  Thus,  in  a  collection  of  papers  with  Isidori, 
Kokotovic.  Kadiyala.  Teel  and  Di  Benedetto,  we  have  laid  out  the  rudiments 
of  a  theory  of  parameter  adaptive  control  for  nonlinear  tracking  and  regula¬ 
tion.  direct  and  indirect  for  systems  with  parameter  uncertainty.  This  result 
has  also  generated  a  large  volume  of  activity  in  the  research  community  and 
has  also  helped  us  understand  adaptive  splining  for  our  CAD  design  package. 
In  future  work  we  will  develop  a  complete  theory  of  Nonlinear  Model  Refer¬ 
ence  .Adaptive  Control.  .Another  area  that  we  will  explore  is  the  question  of 
how  to  identify  nonlinear  models  which  are  presented  in  non-symbolic  form 
using  techniques  from  approximation  theory.  This  will  also  enable  us  to  move 
away  from  exclusively  spbned  approximations  for  our  CAD  package. 

1 .2  Structure  of  Zero  Dynamics  of  Nonlinear  Systems 

Stability  properties  of  zero  dynamics  are  among  the  crucial  input-output 
properties  of  both  linear  and  nonlinear  systems.  Unstable,  or  “non-minimum 
phase",  zero  dynamics  are  a  major  obstacle  to  input-output  linearization 
and  high  gain  designs.  .An  analysis  of  the  effects  of  regular  perturbations 
in  system  equations  on  zero  dynamics  shows  that,  whenever  a  perturbation 
increases  the  system’s  relative  degree,  it  manifests  itself  as  a  singular  pertur¬ 
bation  of  zero  dynamics.  In  this  work,  conditions  are  given  under  which  the 
zero  dynamics  evolve  in  two  time  scales  characteristic  of  a  standard  singular 
perturbation  form  that  allows  a  separate  analysis  of  slow  and  fast  parts  of 
the  zero  dynamics.  The  slow  part  is  shown  to  be  identical  to  the  zero  dy¬ 
namics  of  the  unperturbed  system,  while  the  fast  part,  represented  by  the 
so  called  boundary  layer  system,  desenbes  the  effects  of  perturbations.  It  is 
remarkable  that,  as  the  perturbation  parameter  e  tends  to  zero,  the  bound¬ 
ary  layer  system  becomes  a  linear  system,  whose  stability  is  easy  to  analyze. 
When  this  system  is  unstable  the  perturbed  systems  is  slightly  non-minimum 
phase  and  the  exact  nonlinearity  cancellation  or  a  high  gain  design  should 
be  avoided,  (see  publications  [1]  and  [lOj). 
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1 .3  Approximate  Input-Output  Linearization  by  State 
Feedback 

In  a  collection  of  papers  with  Hauser  and  Kokotovic  we  began  a  study  of  how 
to  enlarge  the  domain  of  applicability  of  nonlinear  control  laws  to  systems 
which  did  not  fit  all  the  assumptions  for  the  rigorous  application  of  the 
theoretical  results.  Our  work  was  strongly  influenced  by  two  important  and 
practical  design  examples:  the  dynamical  model  of  V/STOL  aircraft  and 
that  of  a  ball  and  beam  (modeling  slosh)  in  fuel  tanks  on  aircraft  wings. 
This  work  was  primarily  for  SISO  systems  and  has  been  important  in  that  it 
has  spawned  a  large  effort  on  the  part  of  the  research  community  at  large  on 
approximate  linearization.  In  our  own  work  it  has  been  important  in  helping 
us  develop  the  design  CAD  package  In  future  work  we  will  extend  this  with 
researchers  from  the  University  of  Rome  to  MIMO  systems.  The  subtleties  of 
the  theory  in  these  systems  make  this  a  very  challenging  enterprise.  This  will 
be,  to  our  knowledge,  the  first  attempt  also  to  confront  MIMO  robustness 
issues  for  nonlinear  systems  head  on.  One  report  describing  a  robust  version 
of  the  Descusse  Moog  algorithm  for  dynamic  decoupling  is  under  preparation, 
(publication  [ll]  and  one  more  in  preparation). 

1.4  CAD  and  Implementational  Tools  for  Nonlinear 
Systems 

The  chief  drawback  of  the  recent  advances  in  nonlinear  control  has  been  that 
they  have  been  based  on  detailed  analytical  models  of  the  systems  to  be  con¬ 
trolled.  These  analytical  models  are  required,  since  the  design  methodology 
involves  in  a  fundamental  way  differentiation  of  the  functions  describing  the 
dynamics.  The  reason  that  this  has  been  a  drawback  is  that  there  are  in 
practice  a  large  number  of  nonlinear  systems  whose  parts  are  described  by 
tabular  means  or  in  some  instances  from  empirical  observations.  These  prob¬ 
lems  are  especially  acute  in  flight  control,  where  the  aerodynamic  or  wind 
tunnel  data  is  available  only  at  discrete  points  in  the  flight  envelop  or  in 
the  instance  of  Are  control  for  helicopters  where  a  large  number  of  the  non¬ 
linear  parameters  can  only  be  measured  empirically.  In  collaboration  with 
Integrated  Systems  Incorporated  (ISI)  and  Picatinny  Arsenal  we  have  been 
developing  a  computer  aided  design  package  for  spline  fitting  graphical  data 
and  then  computing  input-output  linearizing  control  laws,  approximate  lin- 
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earization  control  laws  and  also  observers.  The  package  is  written  in  C  and 
involves  adaptation  on  the  order  of  the  spline  fit  as  well  as  the  accuracy  of  the 
approximate  linearization.  The  CAD  tools  being  developed  are  fast  enough 
to  be  real  time  and  the  coefficients  of  the  spline  fit  are  recomputed.  The 
primary  developer  of  the  software  (R.  Kadiyala)  is  also  involved  in  validat¬ 
ing  the  software  on  the  gun  control  models  being  developed  by  ISI  and  the 
fire  control  group  of  Dr.  Coleman  at  Picatinny  .Arsenal.  We  are  also  actively 
exchanging  ideas  and  software  with  a  group  at  University  of  California,  Davis 
under  Prof.  .Arthur  Krener  for  approximate  linearization. 

Thus,  one  of  the  mo.st  important  goals  on  this  grant,  also  continuing 
forward  with  the  next  grant  is  to  develop  at  least  at  a  conceptual  level  user 
friendly  tools  for  tionlinear  control,  which  contain  on  the  one  hand  recent 
advances  in  the  theory,  but  on  the  other  hand  also  take  advantage  of  recent 
advances  in  workstations  to  provide  graphical  and  symbolic  visualization  of 
simulations.  Our  software  has  incorporated  graphical  depiction  of  our  control 
laws  on  Sun  workstations.  This,  we  believe,  is  essential  to  allow  for  rapid 
prototyping  of  new  nonlinear  and  adaptive  control  laws.  The  systematic 
development  of  the  software  in  C  with  a  good  user  interface  are  current  topics 
of  research.  What  has  begun  as  an  off-line  C.AD  tool  design  effort  has,  owing 
to  the  development  of  computer  hardware,  become  an  attractive  option  for 
real  time  control:  consequently  the  real  time  aspects  of  the  computations  are 
our  future  priorities,  (publications  [.5]  and  [13]) 

1.5  Robust  and  Adaptive  Nonlinear  Output  Regula¬ 
tion 

A  new  topic  of  excitement  in  the  nonlinear  control  design  literature  has  been 
the  development  of  techniques  of  out  pul  regulation  for  nonlinear  systems.  In 
contrast  to  the  work  on  input-output  linearization  by  state  feedback  which  is 
the  nonlinear  analog  of  a  zero  cancelling  control  law.  these  methods  do  not 
need  the  underlying  system  to  be  non-ininimum  phase.  Our  research  on  this 
grant  has  been  aimed  at  understanding  the  robustness  of  these  control  laws 
to  parametric  uncertainty.  In  the  instance  that  the  parametric  uncertainty 
is  too  large,  an  adaptive  scheme  is  proposed  with  slowly  varying  parameter 
update  to  achieve  asymptotic  regulation. 

Related  research  in  this  area  concerns  enlarging  the  domains  of  attraction 
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for  these  control  laws,  since  they  are  originally  derived  to  be  local  control 
laws.  Pi  key  difficulty  with  the  new  scheme  appears  to  be  very  small  domains 
of  attraction.  In  research  to  date,  we  have  proposed  several  augmented  to 
the  extended  scheme  to  enhance  the  domains  of  attraction,  (publications  [7] 
and  [8]) 

1.6  Sliding  Mode  Control  of  MIMO  Nonlinear  Sys¬ 
tems 

The  problem  of  developing  precise  matching  conditions  for  nonlinear  systems 
which  are  not  linearizable  by  static  state  feedback  has  proved  to  be  a  surpris¬ 
ingly  hard  nut  to  crack.  In  early  work  on  the  grant  ..e  encountered  success 
in  developing  matching  conditions  for  MIMO  systems  linearizable  by  static 
state  feedback.  The  extension  of  these  results  to  either  dynamically  decou- 
plable  MIMO  systems  or  other  more  general  systems  is  not  yet  complete. 

However,  our  earlier  experiments  with  sliding  mode  control  laws  have 
enabled  us  to  understand  solutions  to  stabilization  problems  where  it  may  be 
shown  that  the  underlying  control  system  cannot  be  stabilized  by  continuous, 
state  feedback,  (publication  [14]) 

2  Scientific  Personnel  and  Degrees  awarded 

1.  S.  Behtash  —  Ph.  D.  awarded  .January  1989. 

2.  J.  Hauser  —  Ph.  D.  awarded  .August  1989. 

3.  R.  Kadiyala  —  M.S.  awarded  Docomber  1989.  (Ph.  D.  expected  Dec 

1991) 

4.  A.  Teel  —  M.  S.  awarded  D**n‘ml)«T  1989.  (Ph.  D.  expected  March 

1992) 

5.  A.  K.  Pradeep  —  Ph.  D.  exp<‘<.tc<l  December  1991. 

6.  Prof.  P.  V.  Kokotovic  —  Visiting  Professor,  Fall  1988. 

7.  Prof.  M.  D.  Di  Benedetto  —  Visiting  Professor,  Fall  1990. 

8.  Prof.;  S.,  S.  Sastry 
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Abstract 

In  tkif  dusartation,  wt  tmbark  on  a  projtct  to  makt  recant  tkaoivtical  advances 
in  geometric  nonlinear  control  into  a  praetieabk  control  dooign  methodology. 

The  method  of  inpnt*ontpnt  Uneariration  by  atate  feedback  provide!  a  natural 
framework  to  deaign  contrdlerf  for  lyitema,  rack  aa  aircraft,  where  ontpnt  traddng  rather 
than  stabilization  is  the  control  objective.  Central  notions  inclnde  relative  degree  and  zero 
dynamics.  Ronghly  speaking,  the  relative  degree  of  a  system  is  the  dimension  of  the  part 
of  the  system  that  can  be  inpnt*otttpttt  linearized  and  the  zero  dynamics  are  the  remaijung 
(nnobservable)  dynamia.  Systems  with  exponentially  stable  zero  dynamics  are  anaiogotu 
to  minimum  phase  linear  systems  and  can  be  controlled  to  track  a  rich  class  of  output 
trajectories  with  internal  stability. 

While  investigatinf  the  nse  of  these  methods  in  the  control  of  the  V/STOL  Harrier 
airaaft,  we  noticed  that  the  small  forces  produced  when  generating  body  moments  caused 
the  airaaft  to  have  an  unstable  zero  dynamics,  i.e.,  to  be  nonminimum  phase.  However, 
if  this  coupling  were  zero,  then  the  aircraft  could  be  iaput<outpu.'  linearized  with  no  zero 
dynamics.  In  other  words,  a  small  change  in  a  parameter  resulted  in  a  significant  change  in 
the  system  structure! 

With  this  observation  as  the  driving  force,  this  dissertation  studies  the  effects  of 
system  perturbations  on  the  structure  of  the  system  and  develops  methods  for  tracking 
contrdla  design  based  on  approximate  systems. 

After  reviewing  the  basics  of  geometric  nonlinear  control,  we  show  that  small  reg¬ 
ular  perturbations  in  the  system  can  result  in  singular  perturbations  in  the  zero  dynamics. 
We  give  asymptotic  formulas  for  the  resulting  fast  dynamics. 


i 


•  • 
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Next,  we  devdop  teduuqees  for  tieddsg  control  design  fm  systems  that  do  not 
have  a  well  defined  relative  degree.  Using  an  approximate  system  with  a  well  defined 
relative  degree,  we  design  tracking  controllers  that  gnarantee  approximate  tracking  for  the 
true  system.  This  approach  is  shown  to  he  superior  to  the  usual  Jacobian  linearization 
method  on  a  simple  baU  and  beam  system. 

Returning  to  the  airaaft  control  problem,  we  use  a  highly  simplified  planar  VTOL 
airaaft  model  to  illustrate  the  (slight)  nonminimnm  phase  characteristic  of  these  systems 
and  develop  a  controller  to  guarantee  approximate  tracking.  We  also  develop  a  formal 
theory  for  this  class  of  systems. 


Toward  Larger  Domains  of  Attraction  for  Local  Nonlinear  Control 

Schemes  * 

Andrew  R.  Teel 

Department  of  Electrical  Engineering 
and  Computer  Sciences 
University  of  California 
Berkeley,  CA  94720 


Abstract 

This  paper  is  motivated  by  the  observation  that  the 
success  of  some  recent  nonlinear  control  approaches  is 
very  sensitive  to  initial  conditions.  The  discussion  in 
this  paper  centers  around  the  recently  developed  non¬ 
linear  output  regulation  theory  of  (l).  The  focus  of  this 
paper  will  be  extending  the  region  of  attraction  for  this 
approach  by  augmenting  the  existing  scheme. 

Keywords.  Nonlinear  Output  Regulation,  Center 
Manifolds,  Domains  of  Attraction. 


1  Introduction 

The  tool  box  for  achieving  tracking  in  nonlinear  sys¬ 
tems  is  growing  rapidly.  One  recent  addition  has  been 
the  nonlinear  regulator  of  (l).  This  solution  allows  the 
control  designer  the  flexibility  of  using  either  output 
feedback  or  state  feedback,  permits  disturbances  to  the 
plant  and  does  not  require  a  well-defined  relative  de¬ 
gree.  However,  it  has  been  found,  through  simulation 
studies,  that  this  control  scheme  can  be  rendered  in¬ 
effective  by  the  choice  of  seemingly  benign  initial  con¬ 
ditions.  In  what  follows,  we  propose  a  means  of  aug¬ 
menting  the  scheme  that  is  provably  convergent  (lo¬ 
cally)  and,  in  simulations,  displays  an  increased  region 
of  effectiveness. 


2  Problem  Statement 

We  begin  by  reviewing  the  problem  statement  and  so¬ 
lution  found  in  (l).  Consider  the  nonlinear  composite 
system 

i  =  /(*)  +  ^(x)u -f  p<r)u; 

w  =  8{w)  (1) 

e  =  h{x)  +  q{w) 

where  *  €  f/  C  R"  is  the  state  of  the  plant,  w  £  W  C 

R'  is  the  state  of  an  (autonomous)  exosystem,  u  €  R'" 
and  e  6  R*.  As  usual,  /  and  the  columns  of  g  and  p 
are  assumed  to  be  smooth  vector  fields  and  h{x)  is  a 

*R«««arch  lupported  in  part  by  th«  Army  under  grant  ARO 
DAAUM-KOSn,  and  NASA  under  pant  NAOS-SO. 


smooth  mapping  on  U.  Also',  s  is  a  smooth  vector  field 
and  q{w)  is  a  smooth  mapping  defined  on  W .  Further, 
it  is  assumed  that  /(O)  =  0,  s(0)  =  0,  h{Q)  =  0,  *(0)  s 
0  so  that,  for  u  s  0,  the  composite  system  (1)  has 
an  equilibrium  state  (x,ip)  =  (0,0)  which  yields  rero 
error. 

We  focus  on  the  following  state  feedback  regulator 
problem:  Given  a  nonlinear  system  of  the  form  (1), 
find,  if  possible,  a  feedback  u  =  a(z,  w)  such  that 

1.  the  equilibrium  z  =  0  of 

*=:/(z)-^ff(r)o(x,0)  (2) 

is  asymptotically  stable  in  the  first  approximation. 

2.  there  exists  a  neighborhood  K  C  x  of  (0,0) 
such  that,  for  each  initial  condition  (x(0),  w(0))  6 
V,  the  solution  of  the  closed  loop  system  satisfies 

Wm  {h{x{t))  +  q(w{t)))  =  0 

f-^OO 


For  the  solution  to  this  problem  in  [l],  the  following 
two  hypotheses  were  made: 

(HI)  the  linear  approximation  of  the  plant  discon¬ 
nected  from  the  exosystem  is  stabilizable.  * 

(H2)  the  point  u;  =  0  is  a  stable  equilibrium  of  the 
exoeystem,  and  there  is  an  open  neighborhood  of  l  < 
point  u;  =  0  in  which  every  point  is  Poisson  stable.  » 
short,  this  assumption  implies  that  the  eigenvalues 
the  linear  approximation  of  the  exosystem  lie  <>**  ‘ 
imaginary  axis.  . 

The  following  solution  to  this  problem  was  then 
veloped: 


Theorem  2.1  (Byrnes,  Isidori)  Under  hypo  * 
(HI)  and  (HI),  the  state  feedback  regulator  proU* 
solvable  if  and  only  if  there  exist  C^(k  >  2)  ^ 

X  =  t(w),  with  w(0)  *  0  end  u  s  e(u)),  vjth  c(0  . 

both  defined  in  s  neighborhood  W*  C  W  ®/®'  ** 
the  conditions 


M>r(w))  ■H(u») 


/(»(«»))+  ,  . 

p(»(w))c(w)+p(»H'“^ 

0 
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1.  The  proof  relies  on  center  manifold  theory  and 
constructs  a  state  feedback 

u  =  a(z,  tt)  =  c(u;)  +  K[x  -  »(ui)]  (4) 

that  is  shown  to  be  a  solution  of  the  state  feed¬ 
back  regulator  problem,  /f  is  a  matrix  of  feed¬ 
back  gains  such  that  the  eigenvalues  of  the  linear 
approximation  of  the  plant  (disconnected  from  the 
exosystem)  have  negative  real  part.  The  manifold 
z  s  )r(u;)  is  seen  to  be  an  error-zeroing  manifold 
that  is  rendered  invariant  by  the  control  u  s  e(w). 
Solving  the  state  feedback  regulator  problem  re¬ 
duces  to  solving  for  the  the  mappings  t  =  t{w) 
and  u  =  c(u»). 

2.  A  very  useful  observation  was  made  in  [2]  that 
these  mapping  could  be  approximated  up  to  ar¬ 
bitrary  order  and  still  achieve  approximate  track¬ 
ing.  This  result  also  followed  from  center  manifold 
theory.,  (The  same  observation  was  made  in  [4]  re¬ 
garding  a  similar  solution  to  this  same  problem.) 
This  observation  makes  actual  application  of  the 
nonlinear  regulator  theory  more  feasible. 

>  Augmenting  the  Solution 

thi.<  section,  we  propose  to  augment  the  solution  to 
«tate  feedback  regulator  problem  given  in  (l}.  We 
.-in  by  motivating  this  augmentation  with  an  appli- 
i<<n  example  initially  studied  in  [3]  and  later  in  (2], 
And  (S). 

unsider  the  well-known  ball  and  beam  example. 

'  It  dynamic  system  can  be  modeled  by 

*1  = 

zi  =  xtzl  -  C7sin(zz) 

*3  S  z<  (5) 

Z4  S  U 

y  =  h(z)  =  *i 

"here  Z|  is  ball  position,  zj  is  ball  velocity,  za  is  the 
^ngle  of  the  beam,  and  z^  is  the  beam’s  angular  veloc- 
|>.v  (For  a  derivation  of  these  equations,  see  (3).  For  ail 
the  simulation  results  mentioned  in  this  paper,  the 
'^^celeration  due  to  gravity  was  taken  to  be  the  nor- 
oulized  value  1.  The  magnitude  of  the  signal  to  track 
then  relative  to  this  normalized  value.) 

The  task  at  hand  is  to  cause  the  ball  position  Z| 
(at  least  almost)  track  a  sinusoid  produced  by  the 
'’’'osysiem 

u»i  =  -Auij 

til.  =  Au/i  (6) 

,(u/)  =  -uii 

presented  in  [2],  (4)  and  [5],  approximating  the 
'"■'nifold  to  either  first  or  third  order  yielded  nice  ap- 
’  '"'’'imaie  tracking  results.  However,  as  discovered  in 
the  region  of  attraction  to  this  manifold  could  be 


very  small.  Further,  the  region  of  attraction  did  not 
seem  to  improve  as  the  order  of  the  approximation  was 
increased.  The  simulation  results  of  section  4  ^lell  this 
out  in  more  detail. 

With  this  in  mind,  we  propose  augmenting  the  so¬ 
lution  to  the  state  feedback  regulator  problem  in  the 
following  manner.  First  consider  the  composite  system 
in  expanded  form 

z  «  Az-b  £lu-P  Pto  +  d(*,w,u) 

v>  =  SvB  +  ’ 

We  retain  hypothesis  (H2)  and  modify  hypothesis  (HI) 
in  the  following  way: 

(Rla)  the  pair  (A,B)  is  controllable. 

Now  augment  the  exoeystem  with  the  following  dy¬ 
namics 

y  as  Ay-f  Bff,y  +  flAf(c)y  ... 

€  sa  0 

where  y  €  V  C  R",  <  6  iP  C  R'"",  and  Af(t)  is  an 
m  X  n  matrix  with  smooth  entries.  We  subject  the 
augmented  exosystem  to  the  following  hypotheses; 

(H3)  is  such  that  all  of  the  eigenvalues  of  (A  + 
BK,)  have  zero  real  part. 

(R4)  ((0)  is  chosen  sufficiently  small  and  such  that 
all  of  the  eigenvalues  of  (A  4-  BK0  -f  BA/(((0)))  have 
negative  real  part. 

The  initial  conditions  for  y  will  be  specified  in  the 
following  theorem  which  is  analogous  to  the  theorem 

ofll). 

Tbnomm  3.1  ffader  hypofkeiei  (Hit),  (Ht),  (B5), 
(H4),  the  title  feedback  rt§tlttor  problem  it  tolvable  if 
there  entt  C^{k  >  2)  mappiayi  z  s  d>(v>,V,f),  vith 
0(0, 0,0)  a:  0  end  u  as  d(w,y,f),  with  d(0,0,0)  as  0 
both  defined  in  a  neighborhood  W*  xY*  x  E*  CW  x 
Y  X  E  of  OxQxO,  tttitfping  the  eonditiont 

^f(w)  +  |J(Ay  +  BK,y  +  BM(f)y)  as 

/(0(u».  y.  <))  +  y(0(w,  u<  <))«<(«'.  V,  <)  +  p(0(«',  y,  <))«; 

*(0(w.  y.  0)  +  ?(»)  -  ^(y)  =  0 

(0) 

end  y(0)  w  tueh  thtt  |z(0)  -  0(ui(O),y(O),t(O))|  it  ttfi 
fieiently  tmtli 

Proof  The  proof  follows  the  proof  of  theorem  2.1. 
Accordingly,  assume  the  conditions  (9)  are  satisfied 
and  consider  u  a  possible  solution  the  state  feedback 

o(z,  w,  y,  f )  =  d(w,  y,  <)  -f  /f |z  -  0(u/,  y,  c)l 

where  all  the  eigenvalues  of  (A  4-  BK)  have  negative 
real  part.  The  existence  of  t  K  such  that  this  is  true 
follows  from  hypothesis  (Hla).  We  now  check  that  this 
state  feedback  is  a  solution  to  the  state  feedback  reg¬ 
ulator  problem.  Requirement  (i)  is  satisfied  because 
o(z.  0,0,0)  s  Hz.  From  hypotheses  (H2)  and  (H3), 
the  overall  composite  system  can  be  transformed  into 


coordinates  for  which  center  manifold  theory  directly 
applies.  Since  a(rl>{w,  y,  f ),  w,  y,  t )  =  d{w,  y,  <),  by  con¬ 
struction  X  =  t^(w,  y,  c)  is  such  a  manifold  in  the  orig¬ 
inal  coordinates.  Also,  by  (9),  the  error  is  given  by 

e(0  =  M*(0)  -  V.O)  +  MlKO) 

It  follows  from  hypotheses  (H2),  (H4),  the  choice 
of  K  and  the  triangular  structure  of  the  compos¬ 
ite  system  that  the  point  (z.ui.y.r)  =  (0,0, 0,0)  is 
a  stable  equilibrium  for  the  composite  system.  So, 
for  sufficiently  small  (z(0),  w(0),  y(0),  r(0)),  the  so¬ 
lution  (x{t),w{t),y{t),i(t))  remains  in  an  arbitrar¬ 
ily  small  neighborhood  of  (0,0, 0,0)  for  all  (  >  0. 
(Notice  .r(0)  and  u;(0)  sufficiently  small  are  provi¬ 
sions  from  the  problem  statement,  <(0)  sufficiently 
small  follows  from  hypothesis  (H3),  and  y(0)  suffi¬ 
cient  I  v  small  follows  from  the  choice  of  y(0)  such  that 
|/(0)  -  t(u(0),y(0),t(0))|  is  sufficiently  small  together 
wiih  /(iM.u'(O),  and  r(0)  sufficiently  small.)  With  this 
sul  iliu  property,  we  can  apply  a  property  of  center 
ii.si.ii  i  Is  yielding  there  exist  real  numbers  Af  >  0  and 
u  •  '’ii'h  that 

•  ■  t(‘*'.y.0(011  <  We"“ll*(0)-i^(«".V.«)(0)ll 

!  t  '  0  Finally,  from  the  continuity  of  h,  together 

o'*  >[  ’thesis  (H4)  and  the  fact  that  h(0)  =  0,  we 
It'-  I'l  e(t)  =  0.  We  conclude  that  this  choice  of 
r-  !’  1  k  solves  the  state  feedback  regulator  problem. 

iti'iiiarks. 

1  III’  manifold  z  s  tl>{w,y,()  is  an  error-teroing 
manifold  in  the  limit  u  y  —  0.  This  manifold  is 
rendered  invariant  by  the  control  u  s  d{w,y,t). 

2  Tlieorclically,  in  terms  of  regions  of  attraction,  we 
do  not  gain  anything  over  the  result  in  [l]  be¬ 
cause  we  are  still  dealing  with  fairly  unspecified 
local  neighborhoods  of  the  origin.  However,  the 
improved  simulation  results  in  some  instances  are 
quite  striking. 

3.  The  reason  for  the  improved  simulation  perfor¬ 
mance  is  that,  with  the  additional  y  states  we  have 
created  an  augmented  error 

e.(0  =  ^(*(0)  +  9(«'(0)  -  MiKO) 

for  which  z  =  4>{w,y,t)  is  an  error-zeroing  man¬ 
ifold.  Dy  quickly  regulating  to  the  manifold  z  =; 
ip{w,y,i),  the  control  steers  the  system  slowly  to 
the  original  error-zeroing  manifold  z  =  t{w).  Reg¬ 
ulating  to  the  manifold  z  =  il>{w,y,()  is  rela¬ 
tively  easy  because  the  system  trajectory  neces¬ 
sarily  starts  close  to  this  manifold. 

4  Examples  and  Simulations 

In  this  section  we  begin  by  presenting  an  example  that 
clearly  demonstrates  the  augmented  solution  proposed 


Figure  1:  Tracking  Results  -  Standard  Nonlinear  Reg¬ 
ulator  for  example  4.2.  (The  dotted  line  (lower)  repre¬ 
sents  the  desired  output.  The  solid  line  represents  the 
actual  output.) 


here.  Then  we  demonstrate  the  usefulness  of  this  aug¬ 
mented  scheme  on  the  ball  and  beam  example.  We  pro¬ 
vide  simulation  results  that  demonstrate  this  Kbeme’s 
ability  to  handle  a  wider  range  of  initial  conditions  in 
the  plant. 

Example  4.1  Consider  the  system 

Z|  =  Zi -f  zj -h  zj -f  Ull 
Zj  =  u 

till  =  -wj  (I®) 

til)  S  till 

e  =  fi(z) -f  ^(ui)  =  zi  -  till 
Augment  the  exosystem  in  the  following  manner; 
ill  =  vi  +  yi 

in  =  -yi  -  yj  -  fiyi  -  <2yz  /JD 

Cl  =  0 

<J  s  0 

with  f  1 ,  (}  sufficiently  small  and  <i  >  tz-  We  now  soltw 
for  the  mapping  z  =  il>{w,y,€)  and  u  =  c(ui,y.t)  ■* 
specified  in  theorem  3.1  and  find: 

=  tiii-fy, 

V’2(«‘'.y.<)  =  -2uii  -  uij  -  (uii -f  yi)* +  y» 
c(w,  y,  f)  =  2ui3  -  ti/i  —  2(ui| -f  yi)(yi  4- W  " 

_y,  _  y,  _  <,yi  -  <2y3 

(*• 

Finally,  we  choose  (yi(0),  yj(0))  in  a  neighborhood 
the  point  (yj,  yj)  given  by 

yj  =  z,(0)-tiii(0)  , 

yj  =  Z7(0)  4- 2un(()) -h  ui3(0) -f  ■ 


Example  4.2  Again  consider  the  ball  and  . 
tern  given  in  (5)  and  the  exosystem  (6)  Augo’ 
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Figure  2;  IVMking  Reaults  •  Augmented  Noniineer 
Regulator  for  example  4.2.  (The  dotted  line  (lower) 
.-('presents  the  original  desired  output.  The  solid  line 


presents  the  actual  output.) 

and 

osystem  in  the  following  manner: 

a(*.  v>,  V,  ()  *  d(u>,  y.  <)  +  /ffx  - 1/>( w,  y,  e)|  (15) 

ill  ^  Vi 

where  K  ntabiiites  the  pair  (A,  B).  First  consider  the 

yj  =  -Cys 

control  (14).  Figure  1  shows  the  inability  of  this  stan¬ 

W  *  y4 

(13) 

dard  solution  to  regulate  to  the  desired  trajectory  from 

ii*  -  tivi  +  <8ya  +  ejys  +  <4y4 

the  initial  conditions 

f.  =  0 

X|(0)  s  15.5 

th  r,  sufficiently  small  and  such  that  the 

equilib- 

x,(0)  w  0 

'iiii  point  y  =  0  is  asymptotically  stable. 

We  cal- 

X3(0)  K  0 

1  lie  a  first  order  approximation  to  the  mappings 

X4(0)  u  0 

-  <^  («i.y><)  and  u«  c(w,y,<) 

>w,y,()  s  uij+yj 

Now  consider  the  control  (15).  Figure  2  shows  the  abil¬ 
ity  of  the  augmented  solution  to  regulate  to  the  desired 

irr.y.f)  a  -Awj  +  yj 
a  ^X'^Wi+lfa 
a  -iA’u/j  +  V4 

. v, ()  a  -^x*wi  +  <i vj  +  <iyj  +  tm  +  uy* 

'il.v,  we  choose  y(0)  in  a  neighborhood  of  the  point 

civen  by 


trajectory  from  the  initial  conditions 


»i(0) 

*a(0) 

*s(0) 

*4(0) 


40 

0 

0 

0 


yl  =  x,(0)-u»,(0) 

y;  =  »j(0)  +  Ats,(0) 

Vs  =  »3(0)“  iA’wi(O) 
y;  =  *4(0)+ JA>W,(0) 

'  '"uulation  purposes,  in  the  original  exosystem,  we 
'  =  35,  wi(0)  =  15,  and  «;3(0)  =*  0.  Conse- 

the  task  is  for  the  ball  position,  xi,  to  track 
'*  lisO-  (By  way  of  reminder,  we  continue  to  use 
■  finalued  value  C  =  1.) 

"  "'’w  compare  the  augmented  scheme  (using  first 
^  approximations)  to  the  original  scheme  (using 
'  approximations).  Observe  that,  with  y  5  0, 
'  "f  theorem  3.1  reduce  to  »  and  e  of  theorem 
we  compare  the  control  laws: 

o(x.U')  a  c(«;)  +  Afx  -  x(w)J  (14) 


Note  that  the  initial  error  with  respect  to  the  wiginal 
desired  output  has  been  increased  by  a  factor  of  80. 
Figure  3  shows  the  tracking  error  with  respect  to  the 
augmented  trajectory.  Note  that  the  small  8teady>state 
tracking  error  is  due  to  approximating  the  manifold  to 
first  order. 


5  Conclusion 

In  this  paper,  the  nonlinear  output  regulation  theory 
of  (l)  was  reviewed  and  applied  to  the  ball  and  beam 
example.  It  wm  found  that  for  some  desired  tracking 
signals,  the  region  of  attraction  for  the  error-seroing 
manifold  was  very  small.  An  augmentation  to  the  ex¬ 
isting  scheme  was  proposed  to  handle  a  larger  range 
of litial  conditions.  The  exosystem  was  augmented  in 


such  a  way  that  a  new  manifold  could  be  calculated 
which  paised  arbitrarily  close  to  the  initial  conditions 
of  the  plant  and  asymptotically  decayed  to  the  original 
error-zeroing  manifold.  This  augmented  scheme  was 
demonstrated  in  simulations  using  the  ball  and  beam 
example. 
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Abstract 

The  object  of  this  paper  ia  to  prove  the  stability  of 
an  adaptive  control  Kheme  designed  to  asymptotically 
achieve  output  regulation  for  a  class  of  nonlinear  sys¬ 
tems.  The  solution  proposed  in  [l]  to  the  nonlinear 
output  regulation  problem  is  reviewed  and  the  robust¬ 
ness  of  the  solution  to  parametric  uncertainty  is  ana¬ 
lyzed.  A  standard  adaptive  scheme  is  then  applied  to 
the  problem  and  slowly-varying  results  are  employed 
to  achieve  asymptotic  output  regulation. 

Keywords.  Nonlinear  Outout  Regulation,  Adap¬ 
tive  Control,  Center  Manifold,  Slowly-varying. 

1  Introduction 

The  task  at  hand  is  to  analyte  and  account  for  pa¬ 
rameter  uncertainty  in  the  nonlinear  output  regulation 
problem.  Recent  work  by  Isidori  and  Byrnes  [l]  has 
produced  necessary  and  sufficient  conditions  for  the 
solvability  of  both  the  state  feedback  and  output  feed¬ 
back  regulator  problem  for  a  class  of  nonlinear  systems. 
In  their  work,  the  signals  to  track  are  restricted  to  those 
that  can  be  considered  as  the  output  of  a  Poisson  stable 
exosystem.  Their  analysis  is  based  on  the  local  prop¬ 
erties  of  center  manifolds.  Using  the  work  in  [l]  as  a 
point  of  reference,  this  paper  will  proceed  to  examine 
the  same  problem  in  the  presence  of  parameter  uncer¬ 
tainty.  In  section  2,  we  review  the  nonlinear  regulator 
theory  and  the  solution  developed  in  [l].  In  aection  3 
we  introduce  parametric  uncertainty  to  the  problem.. 
In  section  4  we  lay  the  ground  work  for  our  adaptive 
scheme  by  reviewing  slowly-varying  theory  for  nonlin¬ 
ear  systems.  Finally,  our  adaptive  scheme  is  developed 
in  section  5. 

2  Nonlinear  Regulator  Theory 

The  subsequent  discussion  follows  closely  that  of  (ij. 
riic  class  of  systems  that  will  be  examined  is  those  of 
the  form 

i  =  f(x,6')  +  p{x,e')w  . 

y  =  h(x)  'U 
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where  w  is  the  state  of  an  (autonomous)  exosystem 

w  =  $(w,i*)  (2) 

For  this  system,  we  will  begin  by  considering  €  R'  u 
a  vector  of  knovn  parameters  in  order  to  review  non¬ 
linear  regulator  theory  in  the  absence  of  uncertainty. 
The  control  objective  is  to  have  the  output  track  a  ref¬ 
erence  signal  that  is  the  output  of  the  exosystem  and 
given  by  -9(w(()).  The  plant  (1)  is  auumed  to  have 
m  inputs  and  o  outputs.  The  state  x  of  the  plant  is 
defined  on  a  neighborhood  U  of  the  origin  in  R"..  The 
state  w  of  the  exoeystem  is  defined  on  a  neighborhood 
W  of  the  origin  in  R*.  Further,  /  and  the  columns  of  f 
and  p  are  assumed  to  be  smooth  vector  fields  and  h(x) 
is  a  smooth  mapping  on  {/.  Also,  s  is  a  smooth  vector 
field  and  ^(u;)  is  a  smooth  mapping  defined  on  W,  The 
composite  system  is  then 

X  s  f(x,0*)  +  g(x,0*)u  +  p{x,9*)w 

til  »  (3) 

e  »  A(x)  +  «(u») 

Finally,  it  is  assumed  that  /(O,  •)  =  0,  r(0,  •)  =  0, 
h{0)  =  0,  9(0)  s  0  so  that,  for  u  =  0,  the  compos¬ 
ite  system  (3)  has  an  equilibrium  state  (x,w)  s  (0,0) 
which  yields  zero  error,  independent  of  the  value  of  t*- 
For  the  state  feedback  regulator  problem,  we  seek  a 
state  feedback  of  the  form 

u  =  o(x,w,0*) 

such  that  the  closed  loop  system 

X  =  /(i,0*)  +  p(x,0’)q(x,w,9*)  +  p(x,0‘)«/ 

til  =  3(x,0*) 

e  =  h(x)  +  q(w) 

exhibits  some  stability  property  and  iimi_oo  e(t)  =  0. 
Following  (l),  we  state  the  nonlinear  state  feedback  reg¬ 
ulator  problem  formally. 

State  Feedback  Regulator  Problem.  Given  a 
nonlinear  system  of  the  form  (3),  find,  if  possible,  a 
feedback  u  =  0(2,^, 9’)  such  that 

1.  the  equilibri  im  z  =  0  of 

x  =  /(x.0'}+g(x,&’)o(x,O,0’)  (5) 
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u  asymptotically  stable  in  the  first  approximation. 
i.e. 

+  9{X,  0,  r))|.«o)  c  Cl 

2.  there  exiaU  a  neighborhood  V  C  U  xW  of  (0,0) 
such  that,  for  each  initial  condition  (z(0),ui(0))  € 
V,  the  solution  of  the  closed  loop  system  satisfies 

lim(M*(l))+9(«"(0))*0 

f«s*00 

Under  the  following  two  hypotheses,  statements  con* 
cerning  the  existence  of  a  solution  to  the  state  feedback 
regulator  problem  can  be  formulated; 

(HI)  the  linear  approximation  of  (5)  is  stabilisable. 
(H2)  the  point  w  =  0  is  a  stable  equilibrium  of  the 
exosystem,  and  there  is  an  open  neighborhood  of  the 
point  w  =  0  in  which  every  point  is  Poisson  stable.  In 
short,  this  assumption  implies  that  the  eigenvalues  of 
the  linear  approximation  of  the  exosystem  lie  on  the 
imaginary  axis. 

Byrnes  and  Isidori  state  neceuary  and  sufficient  con¬ 
ditions  for  the  solution  of  the  state  feedback  regulator 
problem. 

Tbeorem  2.1  (Byrnas  and  Isidori) 

Under  kypeMcses  (HI)  sad  (H2),  ikt  slate  fetdiaek 
regulator  preHem  is  salvakle  if  and  onig  if  there  exist 
C*(k  >  2)  mappings  *  ss  ir(w,d*),  with  e(0,tf*)  =  0 
and  u  s  c(w,d*),  with  e(0,d*)  w  0,  hath  defined  in  a 
neighharhoad  W*  CW  of  0,  satisfying  the  eanditians 

ftsiw.n  M  /(ir,d*)  +  p(»,d‘)u» 

g(r,nc(u>,n  (8) 

h(ir(w,f))  +  f(u>)  ss  0 

Remark.  The  proof  relies  on  center  manifold  theory 
and  constructs  a  state  feedback 

u  =  o(x,  w,  f)  =  c(ui,  d*)  +  *•  »(«".  nl  (7) 

that  is  shown  to  be  a  solution  of  the  state  feedback 
regulator  problem.  If*  is  a  matrix  of  feedback  gains 
such  that  the  eigenvalues  of  the  linear  approximation 
of  (5)  have  negative  real  part. 

3  Parametric  Uncertainty 

To  proceed  with  the  dUcussion,  d*  is  now  considered  as 
a  vector  of  anknown  parameters.  Define  the  matrices 

=  [Kl..,,. 

I'hc  following  structural  assumptions  are  now  made 

A  1  The  pair  (A‘,B’)  is  stabiliiaile  for  all  0*  in  a 
’'ftghborkood  of0’. 

A  2  For  all  0*  m  a  ball  around  0’,;  w  =  0  is  a  sta- 
'’Ir  rguilibnum  of  the  ezosyslem,  and  there  is  an  open 
'‘'ighborhood  of  the  point  w  =  0  in  which  every  point 
•'  1‘oisson  stable. 


Before  attempting  to  handle  the  uncertainties  of  the 
plant  and  exosystem  with  adaptation,  the  question  of 
robustness  is  addressed.  In  this  setting,  a  control  is 
formulated  baaed  on  a  model  of  the  composite  system, 
given  by 

z  *  f(x,0*)-ir9{x,0*)u  +  p(x,0*)w 

lit  =  s(w,0*)  (8) 

e  =  A(x)  +  g(w) 

where  d*  is  a )lzed estimate  of  d*.  Assume  the  following; 

A  3  The  estimate  d*  ties  in  ball  around  d* ..  Tkis  ball 
is  aileasi  such  that  the  gains  chosen  to  asymptotically 
stabilise  the  linear  approximation  of  the  model  (discon^ 
nested  from  the  exosystem)  also  asymptotically  stabilise 
the  linear  approximation  of  the  actual  plant  Cdiicon. 
aeefed  from  the  exosystem.) 

A  4  For  (S),  the  conditions  analogous  to  (6)  are  saU 
isfied  by  the  C*  mappings  x  s  w(u),d*),  u  s  c(w,d*), 
for  all  0*  allowed  by  assumption  AS. 

ASA  certainty  eguiualence  feedback  law  of  the  form 

u  s  o(*.tt»,r)  a  e(w,d*)  +  /f[z  -  ir(ui,d*)]  (9) 

is  applied  to  the  actual  composite  system  (S),  where  K* 
is  a  matrix  of  gains  that  stabilises  the  pair  (A*,B*) 
defined  by 

The  stability  of  the  composite  system  with  (9)  as  input 
is  now  examined. 

Theorem  3.1  (Bounded  error  manifold) 

Under  the  assumptions  (AhA5),  the  composite  system 
fS)  with  (9)  as  input  has  a  center  manifold  at  (0,0,  d’), 
the  graph  of  s  mapping 

»=  ♦(u»,d*,d') 

defined  in  s  neighborhood  IV*  x  y  C  IV  x  Ri*  o/(0,d*), 
satisfying  the  condition 

^S(w,0*)  =  /(♦,d*)+p(#,d*)ui 

p(*,d*)(c(ui,d*)  +  A*(#-x)) 

(10) 

As  a  conseguence,  the  solution  (9)  to  the  state  feedback 
regulator  problem  based  on  (8)  yields  bounded  tracking 
error  when  applied  to  fS). 

Sketch  of  Proof.  First  defined  ^  =  d*  -  d*  and  re¬ 
place  0*  by  4  +  d*.  Next  augment  the  exosystem  with 
4  =  0.  From  the  triangular  structure  and  the  assump¬ 
tions  concerning  the  eigenvalues  of  the  plant  discon¬ 
nected  from  the  exosystem,  it  follows  that  the  closed 
loop  composite  system  can  be  transformed  into  coordi¬ 
nates  in  which  center  manifold  theory  directly  applies. 
In  the  original  coordinates,  and  replacing  4  by  d*  -d*, 
this  manifold  is  the  graph  of  mapping  x  =  ♦(u'.d®,  d*) 
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satisfying  the  condition  (10).  (see  {2l  for  details  of  cen¬ 
ter  manifold  theory.) 

Finally,  by  assumption  the  point  (x,  w)  s  (0,0)  is  a 
stable  equilibrium  of  the  closed  loop  composite  system. 
Under  this  condition,  for  sufficiently  small  (c(0),  w(0)), 
bounded  tracking  follows  from  center  manifold  theory 
and  the  continuity  of  h.  □ 

Remarks: 

1.  The  manifold  ^(w,9*,9*)  is  conceptual  and  will 
not  need  to  be  calculated. 

2.  The  preceding  argument  extends  naturally  to  the 
output  feedback  regulator  problem  also  described 
in  (l). 


4  Slowly- Varying  Parameters 

The  question  of  robustnen  is  now  addressed,  under  the 
added  assumption  that  the  parameters  are  allowed  to 
vary  slowly.  Consider  the  state  defined  by  «  s  «  - 
for  fixed  6*  belonging  to  a  compact  set  T. 
The  dynamics  of  the  state  r  are  then 

i  =  $*,$•) 

S  /(x, $•)  +  fit,  0*)a(x,  w,  $*)  +  p(t,  0*)w 

(H) 

From  Theorem  3.1,  for  every  0*  €  f,  the  equilibrium 
point  is0of(U)is  uniformly  uymptotically  stable, 
uniformly  in  the  parameter  0*.  Due  to  this  property 
that  is  uniform  in  the  parameters,  and  the  differentia¬ 
bility  of  9,  the  system  meets  the  requirements  the 
following  useful  lemma  formulated  by  Hoppensteadt  [3] 
and  recently  restated  by  Khalil, Kokotovic  Id]. 


Lenama  4.1  (Hoppensteadt)  Then  tsitiM  s  £ys- 
psnov /saefiea  W(z,d*,d*)  sack  (ket 

«*i((*i)  <  r)  <  «j(W) 
W,(t,0*,t)F(t  +  *,*,0»,0*)  <  -esdri) 
F,.(r.d*,tf’)|  <  c, 

<  c* 

(12) 

/or  all  r  €  Rr  =  {r  €  R"  :  |r|  <  r)  end  {0*,0*)  € 
r  X  r,  whtre  «i(.),  X2(.),  K3(.)  art  stnetly  inertastng 
fanetiont  sad  C|  end  cj  art  nonntgaUvt  constants. 


With  this  Lyapunov  function  in  hand,  the  slowly- 
varying  analysis  proceeds  in  the  following  way..  Allow 
0*  to  vary.  The  dynamics  of  the  state  r  are  now 


f(i-l-*,9,0\0’)-f^0* 


(13) 


Consider  now  the  Lyapunov  function  of  Lemma  4  1  and 
take  Its  derivative  along  the  trajectories  of  the  system 


(13). 

W  =  W,i  +  W,> 

=S  W,f(x -I- ♦,♦,«•,«•)  + (W,. -U'  •. 

w  <  -»C3(l*|)  +  d,|j‘| 

<  -x(W)-».d,|j*| 

where  *  as  nj  •  and  di  ss  c,  +  C3*upr( 
show  that  s  is  stable  for  small  |r(lo)|  and  sii'ih' .. 
small  |fi*|  observe  that  the  set  D  a:  { W  <  «,(,) ;  ^ 
invariant  set  under  the  condition 

l^l<'‘('‘t(«))M 

R  l*(<o)j  <  «r*('‘i(»))  for  »ny  4  <  r,  then,  from  , ; . 
W(ia)  €  D.  Since  D  is  invariant,  (12)  impli(>  i.. 
|x(f)|  <  4 ,  Vf  >  <0.  In  addition,  if  d*,  —  0  as  : 
then  |x(f)|  —  0  ss  <  -*  00  since  ty(s(t),0*(t),ff  i 
as  f  -•  00.  Note  that  we  can  incorporate  pl,iiii  , 
rameter  variations  about  0*  into  variations  of  O'  a 
example  of  such  a  variation  might  be  a  slowiy.v.ir>ii,, 
exoeystem. 

Corollary  4.1  Under  tht  assumptions  . 

anffUitntlf  small  initisf  condilions  (z(0),w(0)).  i 
staMitfo/ ikt  composite  system  (S)  under  the  inpii , 
it  roknsi  to  plant  and  exosysfem  paromtitr  oanaih  ■ 
Hat  art  snfficienltp  slow  end  stag  in  s  neiykkerkcx  .i  . 
Ike  neminel  persmeter  vsfse  0* . 

Proof.  Follows  immediately  from  the  previous  Uimi  . 
and  diKUSsion.  □ 

The  previous  discussion  is  now  applied  to  a  genen 
indirect  adaptive  control  scheme.  Consider  the  con. 
posite  adaptive  system, 

X  SI  /(x,0*)  +  flz,0*)o(z,w,0*)+p(sc,0')w 
w  s:  s(u;,d*) 

^  :=  C(t,w,0*,I) 
e  =  h(3c)-t-q(w) 

(15) 

Corollary  4.3  Undtr  tht  assumptions  (AJ’AS),  for 
sufficiemtlf  small  initial  conditions  (x(0),  io(0)),  the 
staMitf  of  Ike  composite  spstem  (S)  nndtr  tht  input 
(9)  is  rohnst  to  parameter  variations  in  the  control  law 
that  are  snfficitntlp  slow  and  stay  in  a  ntifhhorhood 
of  the  nominal  plant  parameter  value  9* .  Namely,  the 
stability  of  (15)  is  achieved  i/supi>i,|(7(x,tv.4*,t)|  is 
sufficiently  small. 

Remark.  Because  the  parameter  update  law  is  a 
function  of  z  and  0*,  some  additional  analysis  will 
be  required  to  guarantee  a  sufficiently  small  bound  on 
•«Pi>i,|C(*.«'.^*.<)l- 

Corollary  4.3  Under  assumptions  (A1-A5),  for  (15), 
for  sufficiently  small  initial  conditions  (z(0),U'(0))  and 
supt>tf\G{x,w,0*  ,i)\  sufficiently  small,  if0*  converges 
tojome  0  then  x  converges  to  ♦(u/,9, 0‘)  and  the  steady 
state  error,  e(i),  of  system  (15)  is  bounded  and  given 
by 

h{9{u.0,e'))  -  h{t{u,e’)) 
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„  This  is  the  cue  of  «•  -*  0  M  <  -  oo  and 
r-  0  so  that  i(0  —  0  aa  t  -►  oo.  By  definition  of  x, 
^  ^y«ges  to  Then  by  the  continuity  of 

!  <  the  stability  of  the  eompoeite  system,  the  steady 

'  ,„or  is  bounded  and  given  by  h{^{w,  $,$*))  - 
»  .(.r,r)).  0 

Cofoll^y  for  (15), 

Uf  ftfficHnUs  sms//  initial  conditions  (*(0),  to(0))  and 
^,^,„|C(*,u;,d*,f)|s«jflfieientlysinol/,  if  9*  eenwrgca 

limt^mt(l)  -  0 

Proof.  Here  *  converges  to  ♦(w.d’.d*).  Observe  that 
#1  v,  r  ■  ^* )  satisfies  the  same  partial  difierential  equv 
imh  as  v(u;,  d*)  since  't(to,  d*,d*)  is  the  manifold  made 
invariant  by  the  input  (9)  with  d*  =  fi*.  Thus,  from 
iK«  properties  of  center  manifolds,  x  converges  to  the 
v(tt',d*)  of  Theorem  2.1.  Then  Theorem  2.1  implies 
that  /imi.-*o<(0  »  0.  0 

Remark.  Typiially,  it  is  not  possible  to  guarantee 
(orreet  parameter  convergence  a  priori  without  addi* 
lional  assumptions. 

5  Adaptive  Nonlinear  Output 
Regulation 

The  last  result  of  the  previous  section  suggested  that 
if  an  identifier  could  be  constructed  that  guaranteed 
I*  converges  to  t*  then  asymptotic  tracking  would  be 
guaranteed  u  well.  However,  as  is  known  in  the  adap¬ 
tive  literature,  guaranteeing  parameter  convergence  a 
priori  requires  additional  assumptions.  Rather  than 
take  that  approach  here,  a  specific  identifier  will  be 
suggested  that  will  result  in  asymptotic  tracking.  This 
identifier  is  formulated  in  the  mind-set  of  indirect  adap¬ 
tive  control.  Namely,  an  identifier  is  constructed  to 
atimate  plant  parameters  and  then  these  parameters 
arc  used  in  a  certainty  equivalence  control  law.  The 
identifier  uaed  here  is  analogous  to  the  observer-based 
identifier  found  in  [5]. 

Consider  again  the  composite  system  (3)  where  9* 
is  considered  a  constant  but  unknown  parameter  vec¬ 
tor.  The  following  standard  assumption  for  adaptive 
systems  is  made. 

A  8  The  vector  fields  /{x,9')  and  s(w,9')  and  the 
columns  of  g{z,  6’)  and  p(z,9*)  have  the  following  /in- 
esr  parameter  dependence: 

/(*.«•)  =  »:/.(*) 

pA^-n  =  Ef=i 

=  n=i^:  *.(«") 

where  9’,  i  =  1, ...  ,p  ore  unknown  parameters,  which 
appear  linearly,  and  the  smooth  vector  fields  f,(z). 
Si(u’)  are  known. 


Regressors  are  formed  as 

xj(r,u»,u)  s  +  *(*)»*,.,., 

/i»(*)  +  +  Ps.»(*Kl 

xl(w)  =  {si(w),...,s,(u>)J 

where  summation  over  j,  k  is  implied.  Consequently, 
Xg(x,w,u)  €  and  xZi^)  €  R*"'  contain  all  of 
the  nonlinearities  of  the  system.  Now  the  composite 
system  can  be  written  as 

*  = 

ti*  =  xS(w)«*. 

In  what  follows,  the  conventional  notation  for  estimates 
of  unknown  parameters,  9  will  replace  the  previously 
used  9*,.  To  estimate  the  unknown  parameters,  the 
following  identifier  system  is  uaed. 

i  ss  n,(i-r)  +  xj'(*,u>,u)^ 
w  s  w)  +  xl(u')^ 

^  =  -PXs(*,n>,n)Ps(*-*)-PXw(u>)Pw(w-w) 

(16) 

Here  0*  €  R’"'",  Ow  €  R"*'  are  Hurwits  matrices  and 
Pg  6  R"*",  Pw  6  R'**  are  positive  definite  symmetric 
solutions  to  the  Lyapunov  equations 

+  *  “"Innn 

fllPm-^Pm^w  « 

Finally,  p  is  a  small  positive  constant.  Now,  dehnr 
Cf  wi  X  -  X,  Cw  *  li  -  u»,  and  d>  m  d  -  9*.  Then  if..- 
identifier  error  system  becomes 

it  = 

£*  =  flwCw  +  r 

^  =  -PX*(*<n>,u)P,et  -  px*(v))Pu,Ct,. 

Theorem  5.1  Under  the  assumptions  (AI-A6J,  ' 
sufficientlp  small  initial  conditions  (z(0),  u;(0)),,  />•' 
the  composite  system  (9)  snder  (adaptive)  input  /9' 
3p  >  0  of  the  identifier  (16)  such  that 

^  €  Leo, 

9*  ,  Cy,  €  Lqq  n  £*2, 

5.  (r,u/)€  Loo. 

4-  it,iw  €  Leo, 

5.  lim,^oof»(t)  =  limi_ooC*(0  =  0. 

6.  lim,^oo  e(t)  =  h(x(t))  +  q(w{t))  =  0. 

Proof.  Consider  the  Lyapunov  function 

V(er,Cu,,<i>)  =  pe^Ptet  +  pelPe,ee,-¥  (18) 

Taking  the  derivative  of  U  along  the  trajectories  of  (17) 
yields 

1''  =  -pc[€i  -  pc^Cu,  <  0 


268  ECC  91  Eutoptan  Control  Coitfortnet.  Grenoble,  France,  July  2-5  1991 


Hence  0  <  V(t)  <  V{0)  for  all  (  >  0,  ao  that 
y>^,Co,Cm  €  Loo-  Sin<9  V  is  a  poaitive,  monotonically 
decreasing  function,  the  limit  V'(oo)  is  welKdefined  and 

/eo  ,  foo 

V  dr  sspj  {e^e,+eZem)dr  <00 
so  that  e,,em  €  Lj. 

It  is  now  shown  that  p  can  be  chosen  so  that  the 
analysis  of  section  4  holds.  This  will  imply  that  a  re¬ 
mains  bounded.  Consider  the  parameter  update  law 

i  s  ^  s  VI,  u)P,e,  -  pXm{vi)Pwtm 

Sinm  are  smooth,  w  is  bounded  and  € 
Loot  It  follows  that  3m  >  0  and  a  strictly  incrcMing 
function  k\{.)  such  that 

for  all  a  €  Pr  w  {«  €  R"  :  |r|  <  r).  Then  for  the  Lya¬ 
punov  function  of  MCtion  (4),  equation  (14)  becomes 

^  <t  -sj(W)  +  d,l^| 

<  -*3(1*1) +  </iJ>("«  +  «4(|8|)) 

<  -<c(W)  +  pdj«5(W)  +  /jd|m 

where  a  s  kj e and  rc|  s  « icf'.  Now  pick  po 
sufficiently  small  such  that  (a  »  p«dtN|)(.)  is  a  strictly 
increasing  function  of  W,,  Define  k«  s  (a  -  p«dtai). 
Then 

H'<-*e(lV)  +  pd,m 

for  all  p  <  pt .  Now  observe  that  the  set  D  s  { tV  < 
*i(f)}  for  any  f  <  r  is  an  invariant  set  if 

/><  «s(*i(«))/(dim)ap, 

Hence,  if  p  is  chosen  such  that  p  <  min(p«,pi}  then 
D  is  an  invariant  set.  Finally,  if  |x((o)i  <  *i''(*i(9))> 
then  from  (12)  VF((q)  €  D.  Since  D  is  invariant  (12) 
implies  that  |r(t)|  <  «  for  all  (  >  tp. 

Because  x  is  bounded  and  w  is  bounded  by  as¬ 
sumption,  X  is  bounded.  Since  x,  ui  are  bounded, 
X«(x,w,u),  Xu(v>)  are  bounded.  This  implies  e,,  e«, 
are  bounded.  Since  Co,Cr,eo,iiw  €  and  € 
Lp,  limi.^Qo  Co  ^  limi«*oo  *u>  ^  H- 
Finally,  the  convergence  of  the  tracking  error  is 
proved.  Return  to  the  Lyapunov  function  of  (18). 
The  nontrivial  trajectories  corresponding  to  K  =  0  are 
given  by  the  set 

^  ~  {((oifwtd)  '■ 

eg  =  0,f»  =  0,  xl{e,  w,  u)d  =  0,  ^^(u;)^  =  0) 

From  the  definition  of  d,  trajectories  in  this  set  are 
such  that 

xl(x,\t),u)9  =  xr(*.ui,u)e*  . 

xl(u>)9  =  Xw(“')®’ 

From  Theorem  3.1,  <i>(w,9,9")  satisfies  the  condition 
^xl(w)9’  =  xr(*(u/,S.ff*).u-.u)tf‘  (20) 


Further,  by  assumption  (A5),  ir(u;,P)  satisli,^  >, 
ditioB 

s  xTMv>,9),w,u)0 
From  (19),  ir(w,i)  also  satisfies 

|lxX(u/)tf*  =  xl(^(v>,f).u>,u)9‘  , 

Now,  limi..we«  =  lim,..«i«,  =  0  implies 
0.  So  from  Corollary  4.3,  x  converges  to  4(u7i'  •  ■ 
Now  since  w(w,4)  satisfies  the  same  manifold  cqu , 
as  '•(w.fr),  the  properties  of  center  manifolds  n: 
that  X  converges  to  *(w,f).  From  assumption  ( •, 
f(ui)  s  •h(x(vi,i)).  Then,  from  the  continuity  of 

lim«^wc(0  5!  Iimi-«.h(x(0)  +  »(u((t)) 

*  limt.cD  h(*(f ))  ~  h(x{w,  $))  =  i 
O 

6  Conclusion 

This  paper  has  analysed  the  dynamics  of  a  system  h  . 
parameter  uncertainties  in  the  setting  of  nonlinear  r< . 
ulatioB.  For  small  initial  conditions,  the  nonlinear  t< . 
ulator  solutions  were  shown  to  be  robust  to  para 
eter  uncertainties  and  to  slowly-varying  paramcii  ;- 
The  adaptive  nonlinear  regulator  solution  was  cast  im 
this  slowly- varying  framework.  It  wu  shown  then  ih.i 
there  exists  an  identifier  with  sufficiently  small  gam- 
that,  in  conjunction  with  a  certainty  equivalence  con¬ 
trol  law,  yielded  sero  error  tracking  in  the  limit. 
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Abstract 

.’/.I.  ;<d;i«r  discusses  two  schemes  for  the  trucking  eon- 
•  •  •  MIMO  systems  with  parametric  uncertainty  in 
‘.'ynamics.  The  first  approach  is  an  adaptive  ver- 
‘  u  static  feedback  law  for  tracking  control  based 
results  on  asymptoit>.  model  matching  recently 
■  *(i  in  ([DB90b},  [DB90a]).  The  second  approach 
■iJiiptive  version  of  a  (fynamte  preeompensation 
DM87}). 

1  Introduction 

-  III  yetn,  there  hu  been  a  great  deal  of  research 
in  the  adaptive  control  of  nonlinear  systems. 

1  r'  search  has  been  primarily  focused  on  SISO  sys> 
i>  Ml'  and  some  notable  contributions  are  ([KKM89], 
.1KMK891,(S1891.  (PP89]). 

Ill  tins  paper,  we  consider  a  general  MIMO  nonlinear 
P  of  the  form 

x=  /(i)  +  j(x)u  .. 

y=  h(t) 

where  x(()  €  X,  an  open  connected  subset  ofR",  u(()  € 
R"*,  y(0  €  R'.  Further  we  will  assume  that  /  and  the 
columns  of  g,  namely  gt,  are  analytic  vector  fields  on 
R"  and  the  functions  h,  are  real  analytic  functions  on 
R" 

Our  results  on  nonadaptive  tracking  using  static  feed¬ 
back  for  general  MIMO  nonlinear  systems  is  a  by  prod¬ 
uct  of  the  results  of  ([DB90b],  [COB90])  but  have  not 
appeared  in  the  literature,  to  our  knowledge.  The  re¬ 
sults  on  adaptive  asymptotic  tracking  by  dynamic  pre¬ 
compensation  may  be  viewed  as  being  in  the  spirit  of 
([SI89]).  We  also  discuss  schemes  for  adaptive  tracking 
using  static  state  feedback  and  the  general  problem  of 
adaptive  model  reference  control  for  MIMO  nonlinear 
systems.  No  proofs  are  given  in  this  paper;  they  are 
however  available  in  ((08591)). 

"Researdi  supported  in  part  b)  a  McKay  lecturertliip  held  at 
Berkeley  in  September  1990 

'Research  supported  in  part  by  ARO  under  grant  DAAL- 
6Sl\05T2  and  NASA  under  gram  •N.'tG  2.243 
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2  SISO  Adaptive  Linearization 

We  recall  some  reaulU  from  ([SI89])  to  allow  for  a  bet¬ 
ter  understanding  of  the  differences  with  respect  to  the 
MIMO  situation  illustrated  in  Section  3. 

Consider  the  system  (1)  with  one  input  and  one  out¬ 
put.  Let  xo  be  an  equilibrium  point  of  the  undriven 
system,  that  is  /(xq)  a  0  and  that  the  output  is 
zero  at  xo,  i  e.  h(xo)  -  0.  We  will  assume  that  the 
system  (1)  has  s(nc(  relative  degree  y  at  xo  ([Isi89]). 
One  can  then  choose  a  new  set  of  coordinates  given 

by  s  h(x),(,  =  Ljh(x) . =  Lj~^h(x)  and 

i;  €  such  that  dq,  ;  s  0.  In  these  coordinata, 
the  system  (1)  taka  a  ‘normal  form’  which  reveals  the 
sere  dynamics  to  be  of  the  following  form 

ijsqCO.rj)  (2) 

with  r)  6  R'*~''..  The  sero  dynamia  is  said  to  be  expo¬ 
nentially  attractive  to  a  large  ball  in  X  if  the  following 
equation  holds: 

for  l»j|  >  R  (3) 

The  zero  dynamiasatisfia  a  conic  confinuify  conditio* 
in  (  uniformly  in  ij  if 

l7((.»?)-7(0.7)|<t|Cl  (0 

It  may  be  verified  using  a  converse  Lyapunov  argument 
(as  in  ([5189]))  that  asymptotic  tracking  with  bounded 
stata  can  be  obtained  if  the  system  is  exponentially 
attractive  with  conic  continuity. 

The  preceding  rault  has  been  critically  examined  in 
the  literature  but  it  has  not  been  appreciated  that  the 
condition  of  (4)  is  not  a  global  Lipschitz  condition  on 
the  function  q. 

Now,  for  adaptive  tracking,  assume  that  the  vector 
fields  f,g  in  (1)  and  the  function  h  in  (1)  are  un¬ 
known  but  may  be  parameterized  linearly  by  parame¬ 
ters  9'  €  R*  in  the  form 

/(^.«*)=  eLi  <»:/.{*) 

9(T,e')=  Z‘,=i^:9M 

Mr.e')= 


Aiiovvil  luiictions  ot  X.  In  the  equation  (5)  above,  it 
follows  that,  if  some  of  the  0’  are  known,  they  are  re¬ 
placed  by  their  values  The  linearizing  control  laws 
of  the  previous  section  are  replaced  by  their  estimates 
depending  on  the  current  estimate  ${t)  of  6*  in  accor¬ 
dance  with  a  heuristic  known  as  the  certainly  equiva¬ 
lence  prxnctple.  Thus  if  the  "true"  system  is  known  to 
have  relative  degree  y  then  the  control  law  is  given  by 

«  - =4 - +  v)  (8) 

•  L,Ly^h{x) 


.  _ _ _  u<«u  auiitjy  me  conic  continuity  condi¬ 

tions  of  (4).  Also,  assume  that  the  regressor 
has  bounded  derivatives  in  both  its  arguments. 

Then  given  a  bounded  trajectory  yu  vnth  first  y  -  1 
derivatives  alt  bounded  it  follows  that  the  control  law 
of  (6)  with  the  parameter  update  law 

^  =  *  = 

yields  bounded  tracking,  i.e.  y(t)  -»  yM{t)  tvith  all  the 
states  z  bounded,  provided  that  the  state  trajectory  is 
confined  to  X. 


Here  LfLy'h{z),L'fh{x)  stand  for  the  estimates  of 
LjL‘f~'‘h(x),L'jh{x)  derived  by  first  expressing  these 
function  in  terms  of  the  known  vector  fields  fi.gj  and 
known  functions  hk  and  multilinear  products  of  the 
form  6, ..  .0m  and  then  replacing  the  multilinear  prod¬ 
uct  by  an  estimate  of  the  form  6,.  ..0m  ■  We  define  the 
ii'.ultilinear  product  as  a  new  parameter  and  estimate 
II  V  stands  for  the  estimate  of  the  tracking  control  law 
jvcn  by 


<  =  ylf+oniyU^  -  +  ••  ■  +  af,(yw  -  A) 

vnir  that  all  the  L'jhlz)  are  multilinear  functions  of 
macquently,  if  one  defines  0  €  R*  to  be  the  vector 
vil  multilinear  products  of  the  0,  up  to  terms  of 
roe  7,  It  follows  that  the  control  law  of  (6)  is  affine 
'-)  Defining  the  parameter  error  in  ©  to  be  ♦  ;s 
>•)'  and  the  output  error  to  be  <  i=  y  -  y« ,  an  easy 
ilation  yields  that 

e^  +  o,e’-‘+-+a,e5:*^W(z.d)  (7) 

nne  appropriately  defined  W(*,0)  €  R*.  Define 
uiodei)  transfer  function 


Af(s)  = 


1 


s'*'  +  OlSl”*  o, 

'll  augmented  errrr  cj  to  be 


(8) 


■  '  +  (0^(t)A/(s)W(i.0)  -  A/(5)0^(<)W(/,e)) 

(9) 

'■  'inhining  (7)  with  (9)  yields 


e,  =4^A/(s)W(i,0) 


(10) 


Remarks: 

1.,  The  parameter  update  law  is  specified  for  0.  This 
neglects  the  multilinear  dependence  of  terms  in¬ 
side  the  vector.  However,  this  is  necessitated  by 
the  lack  of  a  systematic  theory  of  nonlinear  pa¬ 
rameter  estimation  or  identification. 

2.  Given  the  form  of  the  linear  error  equation,  there 
is  a  large  choice  available  to  us  for  parameter  up¬ 
date  laws.  We  choose  the  normalized  gradient  type 
algorithmof  (11)  here  for  reasons  of  brevity  but  we 
hasten  to  add  that  several  other  normalized  algo¬ 
rithms  (such  as  the  normalized  least  squares)  will 
do  as  well  (see  ([SB89]). 

3  MIMO  Systems 

3.1  IVacking  by  Static  State  Feedback 

Consider  the  plant  P  to  be  a  square,  nonlinear  plant 
of  the  form  (1).  It  is  useful  to  derive  the  tracking  re¬ 
sults  as  a  special  case  of  model  matching  results,  con¬ 
sequently  we  consider  a  model  M  of  the  same  form 
as  the  plant  with  state  z  €  Xm  open  C  R"**  and 
with  f,g,h  replaced  by  fM,gM,hsi.  We  will  need 
to  assume  that  the  f,g,fu,gt4  are  analytic  vector 
fields  and  that  h,hM  are  analytic  functions.  Define 
C(x)  :=  span{yi(r), . .  .,ym(*)}  (over  the  ring  of  ana¬ 
lytic  functions)  and  assume  that  the  dimension  of  the 
distribution  G  is  m  for  all  z^  X.  The  notation  ysrCO 
is  used  to  mean  the  output  of  the  model  starting  from 
state  zo  ot  0  if  there  is  no  need  to  highlight  the  depen¬ 
dence  on  the  initial  state. 

An  extended  system  is  associated  with  the  plant 
and  model  as  follows: 


■uupnient  to  denote  the  filtered  regres.sor 

'  I  a  0)  by 

Wi(i.0)  :=  M(s)W(x,0) 

2.1  Adaptive  IVacking 
'  the  system  of  (1)  with  the  vector  fields  f,g 
fii’iction  h  parameterized  as  in  (5).  Assume 
'//••Icm  can  be  globally  converted  into  normal 
’‘liiinles  on  A'.  Further  assume  that  the  zero 
f  the  system  are  exponentially  altractne  in 


x®=  /®(*®)  +  y(x®)u-fp(x^)u 
y^=  h^(z^) 


with  state  (x®)^  (i^,r^)  €  X  x  Xm,  inputs  u.v 

and 


/V)  = 


/(*) 

.  /«(-) 


j(«) 

0  ' 


,h^(x^)  =  h(x)  -  /iA/(r) 
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Further  define 

Also,  define  the  dynamical  system  with  state  i®,  input 
u  and  output  y®  described  by  the  triple  (/®,  to 
be  t.  Now,  consider  a  point  x§  =  (xo,  ^o)  which  is  an 
equilibrium  point  of  and  also  produces  zero  output 
for  the  system  E,  i  e. 

/®(xf)  =  0  h®(x®)  =  0 
Now,  assume  that 

Assumption  A1  :  (Regularity  of  E) 

Xq  is  a  regular  point  for  the  zero  dynamics  algorithm 
applied  to  E  (regular  in  the  sense  of  ([Ui89],  page  304). 
Actually,  the  assumption  A1  is  a  sufficient  condition 
in  order  to  apply  the  zero  dynamics  algorithm  to  the 
system  E  around  x^.  Let  Mk  denote  the  submanifold 
defined  at  step  k  of  the  algorithm  and  M'  denote  the 
zero  dynamics  manifold  obtained  at  th<.  conclusion  of 
the  algorithm;  there  further  exists  a  unique  smooth 
control  uo  ■  Aif*  -*  R"*  so  as  to  make  M*  invariant, 
i.e.  /^(x^)  +  §{x^)uo  is  tangent  to  A/*.  The  vector 
field  /^(x®)  +  j(x®)uo  restricted  to  M'  is  referred  to 
as  the  zero  dynamics  of  E.  It  can  also  be  shown  that 
M*  can  be  expressed  in  a  neighborhood  of  xf  m 

AY*  =  (x®  €  X  X  Xm  :  H'(x^)  =  0} 

for  some  function  H*,  The  following  theorem  uses 
the  procedure  of  the  zero  dynamics  algorithm  to  solve 
the  model  matching  problem  as  follows  ([DB90a], 
(CDB90]); 

Theorem  3.1  Stable  Model  Matching 
Consider  the  system  of  (IS)  and  assume  (hat  there  ex¬ 
ists  an  x§  such  that 

1.,  A1  holds, 

S.  E  IS  minimum  phase  at  xf ,  and 

S.  span  {^x^)}  C  T,eMs+  span  {j(x^)}  in  a 
neighborhood  of  x§  in  AY*,  for  all  k  >  0. 

Then,  there  exist  neighborhoods  U  of  xq  and  Um  of  xq 
and  an  integer  v  such  that  the  compensator  Q  defined 
by 

X=  a{x,x)  +  b(x,x)v 
u=  c(x.x)  +  d(x.x)v  ' 

for  appropriately  defined  analytic  a,b,c,d  and  x  €  R*' 
and  a  function  F  :  U  x  Um  — *  R*'  and  L  €  R+  such 
that 

a)  If  v{t)  =  0  <Aen  the  point  (xo,Xo  ■=  A'(xo.ro)) 
an  asymptotically  stable  equilibrium  point  of  the  closed 
loop  PoQ,  I.e.  of  the  system 

x=  /(i)  +  9(x)c(x,x) 

X=  o(x.x) 

Y/|v(<)|  <  L  for  allt>0  then 

hm  y^®‘^(x,  F(x,  e),t)-  ysiieJ)  =  0 
!-•<» 

for  all  (i, :)  X  Um 


Remarks: 

1,  In  view  of  the  propositions  of  ((Isifi'j',  \ 
pendix  B.2)  the  fulfillment  of  (a)  above  kh.,,  , 
tees  that,  given  <  >  0,  there  exist  6,l\  s,. 
that  if  |(x(0),x(0))|  <  6  and  |u(.)|  <  a  ,  ,1 

|(x((),x(f))l<«forallt>0. 

2.  The  proof  of  the  preceding  theorem  in  ilic 
mentioned  references  is  constructive  and  iIk.'  .... 
pensator  may  be  shown  to  be  of  the  form 

X=  fM{x)+9M{x)'> 
us=  u(x,x,u) 

initialized  at  xo  =  *0.  Xo  =  F(xo,  zq)  =  zu  a 
a  consequence,  we  have  that  x(0  =  s{t)  and  on 
may  define  the  control  law  in  terms  of  x^  alun. 
rather  than  x,  x.  <  as 

u(x®,u)  :=  u*(x®,i;)  +  M'\x^)KH-{x^)  (l^, 

where  AY(x®)  €  R’"’*'"  :=  dH-(x^)g{x^)  aiui 
u*(x®,  v)  uo(*^)  +  ui{i®)u  is  the  unique  solu¬ 
tion  u  of  the  equation 

dH'(x^)U^{x^)  +  j(x^)u  +  p(x^)u)  =  0 

so  that 

«o(x®)  =  M-\t^)d}r{x^)f^{x^) 
and 

u.(x«)  =  AY-‘(x®)dH*(x^)p{x^) 

Further  K  €  fi"'*'"’  is  chosen  to  stabilize  part  of 
the  system  dynamics  as  specified  below. 

Let  Xq  (respectively  xf )  be  an  equilibrium  point  of 
P  (respectively  E)  such  that  h{xo)  =  0  (respectively 
/i(xf)  =  0).  Then,  the  following  assumptions  aie 
equivalent. 

Assumption  Al'  ;  (Strong  Regularity  of  E) 

E  IS  right-invertible  and  (Xq  =  0)  is  a  strongly 
regular  pair  for  E  (strongly  regular  in  the  sense  of 
(DBG90]). 

Assumption  A2  :  (Strong  Regularity  of  P) 

P  is  right-invertible  and  (xo,  y  =  0)  is  a  strongly  regulai 
pair  for  P. 

Weaker  assumptions  than  A2,  e  g  the  regularity  hy¬ 
pothesis  of  [Isi89]),  p.  302  are  also  sufficient  for  our 
purposes.  The  implications  of  these  assumptions  are 
discussed  in  ((DBS91]).  It  is  shown  in  ([DB90b])  that 
if  the  hypothesis  Al  of  Theorem  3.1  is  replaced  by 
the  assumption  A2  above,  then  one  can  construct  a 
local  change  of  coordinates  (^,17,  x')  =  ♦(x.x)  with 
z‘  ■=  Z  -  Zo,  ♦(xo,  zo)  =  0  and  (  =  H*{x,  z)  such  that 
the  plant  with  the  controller  of  equation  (14)  has  the 
form 

--  /m(x' +  Xo)  +  9m(x' +  xo)<-' 

1]  = 


(1C) 
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The  dependence  of  the  matrices  W2,  IVs  on  the  data, 
and  for  that  matter  on  W\  above,  is  involved.  The 
equations  (24)  are  afhne  in  4  as  a  consequence  of  the 
linear  parameterization  of  the  control  law  by  the  un¬ 
known  parameter  6. 

We  are  not  as  yet  able  to  give  a  stability  proof  for  a 
parameter  update  law  derived  on  the  basis  of  a  com¬ 
posite  Lyapunov  function  involving  the  system  of  equa¬ 
tion  (24)  and  an  equation  for  However,  there  is  one 
important  special  case  for  which  an  adaptive  scheme 
can  be  derived  and  this  is  the  ease  when  the  function 
=  0.  It  can  be  shown  from  Theorem  4.2  on  page  272 
of  (I3i89]  that  this  assumption  is  satisfied  if  the  plant 
P  can  be  linearized  by  static  state  feedback.  This  con¬ 
dition  is  slightly  weaker  than  the  condition  that  P  has 
vector  relative  degree  (in  which  case  it  can  be  both 
linearized  and  decoupled  by  static  slate  feedback) 

Theorem  3.3  Static  State  Feedback  Adaptive 
Tracking 

Consider  the  system  of  (U)  end  the  model  of  (19) 
With  ike  assumption  tkai  9)((,q,z')  s  0  Assume 
that  At  stove  holds  and  tAat  P  is  etponentially 
attroctive  on  X.  Also  assume  that  ike  vtcior 
fields  i('(^.  q,  3').  v,6)  are  Lips. 

cAitr  continuous  in  ihttr  variables  on  X  n  Xm  x  R"  x 
R*'.  Further,  assume  that  Wj  has  hounded  derivatives 
with  respect  to  their  arguments. 

Then,  under  reparomelenxation  unth  the  control  law  of 
(S3),  assume  that  the  system  can  he  expressed  as  in 
equation  (Sf). 

Then  there  easts  a  choice  of  parameter  update 
law  for  4  such  that  the  control  law  of  (S3) 
yields  asymptotic  tracking,  with  hounded  states  when 
suPi>o(lyM(OI,  -M|yw ‘(01)  <  61  <tnd  |*ol.l*(0)|  < 
Sj  and  the  trajectory  of  e  X  x  X^.  The  proof  of 
the  theorem  constructs  the  parameter  update  law. 


with  at,  fit  analytic  functions  of  i  sud,  i 
coupling  matrix  of  E*  with  the  control  |„u  ,  > 
the  form 


A,(x)  = 


0 

hf(x)  0 


Step  3  There  exist  qt  columns  of  Ai(x)  (v. 
of  generality  the  first  qt)  with  two  or  more  i,.., 
ements  Put  an  integrator  in  series  with  (/*  on. 
mg  input  channels,  i.e.  define  the  dynannr . ,,  , 
of  Et  composed  with  (25)  as 


Cl  =  V| 


for  1  =  J,  .  .,qt.  Let  E*+i  be  the  new  system  ol,; 
by  composing  E*  with  (25)  and  (26)  and  return  1 
1  with  k  —  IcA-l  and  »•  *-  (x*)  U{Ci) 

If  the  original  system  is  right-invertible,  then  ilu  ), 
dure  converges  in  a  finite  number  of  steps  to  .1  .s>s;, 
denuted  E*,  having  vector  relative  degree  (rf, 
U’.-9*,h*)  be  the  triple  characterizing  E',  j 
(x,<;)  its  state,  u*  its  input  and  y*  its  output  ( 
struct  a  local  change  of  coordinates  ^(x*)  =  {{,  , 

C  =  coi(Ci)  by  setting 


col  (h, •(*•), L/.h, *(*'), 

-  col  {CuVs . Vr>) 

and  using  some  complementary  coordinates  ij  1  lu 
E*  takes  the  standard  form  ([IsiSS),  pg.  240): 


’)=  9(C>q)  +  p(C,q)u‘ 

{1=  Vs 


Vr> 

Cr  J  =  AJ  (C ,  q)  +  Xl'Jl  T/)u; 
y.‘=  Vi 


(27) 


3.3  Adaptive  Ttacking  by  Dynamic 
State  Feedback 

We  now  turn  our  attention  to  tracking  by  dynamic 
state  feedback.  Several  algorithms  have  been  proposed 
in  the  literature  for  this  problem  and  we  now  recall  the 
one  of  ({DM87]).  We  change  notation  slightly  to  refer 
to  the  process  P  as  Eo.  Set  A  =  0  and  x*  s  x. 

Step  1  Let  r,  be  the  relative  degree  of  the  i  th  output 
of  Et.  Define  the  decoupling  matrix  A(x)  to  have  its 
ij  th  entry 

and  denote  its  normal  or  generic  rank  by  st .  If  s*  =  m, 
stop. 

Step  2  If  s*  <  m,  assume  that  the  first  s*  rows  of 
A(x)  are  linearly  independent  at  each  point  of  an  open, 
dense  subset  of  X  (this  can  always  be  achieved  by  a 
permutation  of  the  components  of  the  output).  Apply 
the  static  stale  feedback 

(25) 


for  i  =  1 .....  m  and  of, ((,ij)  =  L,. iJ; " ‘ hf (^- ‘ (f ,  i,) ) 

‘  <  «.3  <  m  and  6f(^,rj)  =  Lj:hf(^-‘({, »;))  loi 
I  $  •  S  At  this  point,  asymptotic  tracking  may  be 
obtained  by  applying  the  following  control  law; 

u‘=(A‘)-‘(-fr»+  ^ 

yui  +“u(yAn  ‘4r«-l)  +  ••■  +  air;(yMl  -  si) 

.  yMm  +  ®n»l(y*fm  ‘  +  •  •  •  +  Omrj;,(yjWm  -  s 

where  the  polynomials  s’’^  -f  o,is''?”‘  +  a,,*  are 

all  Hurwitz 

Prior  Information  for  Adaptive  Control 
One  assumes  that  the  true  system  is  right  invertible 
Now,  the  variables  are  all  functions  of  the 

unknown  parar^ever  8* .  The  functions  ak,0s  are  func¬ 
tions  of  8*.  To  estimate  these  one  needs  knowledge 
of  the  relative  degrees  of  the  system  E*  at  every  step 
in  the  procedure  above  and  in  particular,  the  vector 
relative  degree  of  the  system  E*  namely  (r{, . .  ,r^) 
Also,,  the  integers  St.gt  tepresenline  thf  rank  of 


U  =  Ot(x)  +  0t{x)v 
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hero  ihc  matrix  A  is  rendered  HurwiU  by  appropriate 
rhoicc  of  l<  in  (15).  The  states  (  contain  in  particu- 
\  ,1,0  output  errors  as  some  of  their  entries.  Also 

,hr  functions  qi  and  p,  satisfy  some  extra  conditions, 

namely 

=  0  ^(0.0.0)  =  0 


Pi(0,»?,  r')sO 

Further,  the  dynamical  system 

s'  =  /m{s'  +  So) 

ij= 


07) 


7-  =  0  closed  loop  system  PoQ  is  asymp¬ 

totically  stable  with  equtlibnum  point  zq. 

2.  When  8up,>o(lyi»r(0l.  -.|y»r‘(0l)  <  and 
|zo|  <  62  then 

lim  =  yuit) 

•  •^00 

3.2  Adaptive  Static  State  Feedback 
Dracking 

In  this  section  we  consider  models  of  the  form  of  equa¬ 
tion  (1)  with  the  added  feature  that  the  dynamics 
of  the  plant  depend  on  certain  unknown  parameters 
€R',  i.e. 


represents  the  zero  dynamics  of  t  and  the  system 


z=  nx,0‘)  +  g{z,e-)u 
y=  h{z,8-) 


V-t(0,q,Q)  (18) 

represents  the  zero  dynamics  of  P.  The  zero  dynamics 
manifold  of  t  is  now  given  by 

A5r-  =  {(<,r,,r')|f  =  0) 

The  form  (17)  of  the  zero  dynamics  of  the  system  E 
shows  that  it  is  minimum  phase  if  the  zero  dynamics  of 
P  and  the  undriven  model  dynamics  are  asymptotically 
stable.  In  fact,  the  decomposition  (16)  can  be  used 
to  extend  the  proof  of  Theorem  3.1  to  cover  the  case 
where,  instead  of  assuming  the  asymptotic  stability  of 
the  zero  dynamics  of  E,  one  assumes  that  the  variables 
z'  ate  bounded  by  a  sufficiently  small  constant  and  that 
the  zero  dynamics  of  P  is  asymptotically  stable.  This 
can  then  be  usefully  applied  to  solve  trajectory  track¬ 
ing  as  a  special  case  of  the  model  matching  problem 
in  which  the  desired  trajectory  ysr  is  generated  by  a 
model  consisting  of  chains  of  integrators  driven  by  the 
appropriate  derivatives  of  the  yut-.  More  precisely,  de¬ 
fine  p,  to  be  the  essential  order  of  the  i  th  output  of  the 
plant  yi  as  defined  in  ([GM89]).  Then  define  the  model 
to  be  matched  to  have  state  z  s  col(r,,t  =  1,.  ,m) 
with  dynamics 

ill  =  *ij 
i|2  =  *13 

(19) 

*111,  =  u, 
ywi  =  *ii 

The  model  satisfies  the  third  hypothesis  of  Theorem  3.1 
by  the  definition  of  the  pi  and  corresponds  to  = 
V|.  Define  p  =  max,p,. 

Theorem  3.2  MIMO  Asymptotic  Tracking 
Assume  that  A2  above  holds  and  that  P  is  minimum 
phase  at  zo  Then,  there  exist  constants  61,62  and  a 
compensator  Q  of  the  form 

u  =  c(x,yM,yM . y/iJ'‘’)+d(a:,yM,yA/... 

(20) 


The  assumption  A2  of  the  previous  section  is  assumed 
to  hold  for  the  true  value  of  the  plant  parameter.  Car¬ 
rying  forward  the  dependence  on  0  through  the  deriva¬ 
tion  of  the  tracking  control  law  will  yield  the  manifold 
H'(x^,0*)  and  the  control  law  of  (15),  namely 

u(z^,u,tf*)  :=  ,v,0')-i-M-\x^ ,0')KH'(x^ ,0') 

(22) 

The  pnor  information  needed  for  adaptive  control  is 
as  follows;  one  assumes  that  at  each  step  of  the  zero- 
dynamics  algorithm  modified  as  described  above  for 
stable  model  matching,  the  manifold  Mi  6  scribed  as 
the  zero  set  of  the  functions  Hk(x^,0)  satisfies  the  con¬ 
dition  that 

dHi{x^,0)g[x^,0) 

has  a  left  null  space  of  constant  dimension  as  a  function 
of  0.,  This  is  a  sort  of  regularity  hypothesis  on  the  plant 
as  a  function  of  0.  Note  that  the  model  is  assumed  to  be 
known  and  independent  of  0.  Since,  the  parameter  0‘ 
is  assumed  unknown,  we  will  replace  it  by  its  estimate 
at  time  f,  denoted  0(f).  Further,  we  will  assume  that 
the  control  law  can  be  linearly*parameterized  as 

u(x®,v,0-)  =  u(x^.u)-t-  tV,(z^,v,0-)0- 

for  an  appropriately  defined  matrix  lVj(x^,u)  € 
and  parameter  vector  6*  €  R*.  Actually  both  u  and 
Wi  are  affine  in  v.  As  a  consequence  the  adaptive  model 
matching  control  law  is  given  by 

u(x®,u,d(l))  :=  u(x^,u)  -F  W',(x^,u,e)0(f)  (23) 

Denoting  the  parameter  error  ♦  =  .0(f)  -  0*  €  R*  we 
may  use  the  control  law  of  (23)  in  the  system  of  (21) 
to  yield  the  following  modification  of  (16)  (note  the 
non-existence  of  pi  in  particular,  caused  by  the  special 
choice  of  model). 

{=  A(  +  qi{(,ii,z')  +  W2((,Ti,z',v,Q)<^ 

*' =  />/(*'  + Co)  +  yAf(*' +  *0)1* 

t)  =  ili(i.T).  z')  -F  z')v  -F  W3((,t},z',v,Q)4> 

(24j 


such  that 
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Imp  matrices  and  the  number  of  integrators  to 
**'  each  step  in  the  algorithm  described  above 

*  a*»uiucd  known  and  are  independent  of  $.  Also,  the 
of  the  integrators  at  each  step  are  assumed 

I!,  m  litis  information  it  is  possible  to  compute  o*,^* 
J  » function  ofO.  As  in  the  SISO  case,  we  will  assume 

^„enil  linearly  on  a  new  parameterisation  0  of  the 
,.»nowii  parameters. 

th-  adaptive  control  law  follows  by  replacing  0  by 
ft  io  that  by  using  the  certainty  equivalence  control 
A«  a  consequence,  the  normal  form  equation  (27)  are 
modifed  to  have  regressor  vectors  iii)(x*,0)  possibly 
at  every  entry  corresponding  to  the  mismatch  between 
t)  and  @*  Thus  the  error  equations  for  the  tracking 
errors  e,  »  y<  -  VMi  ««  given  by 


#,  » 


wliere 


M{  = 


I'f  +o,n'J”'  + 


M}. 


(28) 


Note  that  all  the  transfer  functions  Afj  are  proper,  sta¬ 
ble  transfer  functions.  We  define  the  augmented  error 
to  be 


e,.  =  e.  -F  (Afi(s)ui\)^(t)  -  Af{{s)(uii6(t)) 

-F  • .  +  {M‘.(s)ui’.)6(t)  -  M;j(s)(ui;.6(f)) 

It  IS  easy  to  see  that  the  augmented  error  is  of  the  form 

Cl.  =  W,{x‘,6)*  (29) 


where 

W'(i'.0)  =  Af}(s)u;',(x*,0)-F  .  -F  A/;_.(j)u;',(x‘,e) 
is  a  filtered  regressor.  Note  the  resemblance  of  this 
error  equation  to  that  in  the  SISO  case. 

Under  the  same  hypothesis  as  in  Theorem  2.1  the  same 
conclusions  hold.  There  is  however  one  difference  in  the 
proof  from  the  SISO  case,  namely  that  the  zero  dynam¬ 
ics  are  indeed  driven  by  the  input  u*  in  the  MIMO  case. 
As  a  consequence,  as  in  the  case  of  the  proof  of  Theo¬ 
rem  3.3  we  need  to  insist  that  the  initial  conditions  of 
the  states  x* ,  the  initial  parameter  error  4(0)  and  the 
tracking  output  j/m  and  their  appropriate  derivatives 
are  small  enough. 
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Abstract 

This  paper  discusses  two  schemes  for  the  adaptive  control  of  classes  of  MIMO 
nonlinear  systems  with  parametric  uncertainty  in  their  dynamics.  First,  the  problem 
of  tracking  a  reference  trajectory  is  considered  and  an  adaptive  version  of  the  input- 
output  decoupling  algorithm  of  [DM87]  for  general  right  invertible  MIMO  systems 
is  proposed.  Then  on  the  basis  of  some  results  of  [DB90a],  [DB90b]  on  asymptotic 
model  matching,  a  scheme  is  presented  for  Model  Reference  Adaptive  Control  and  a 
solution  is  given  for  input-output  linearizable  systems.  Moreover,  the  non-adaptive 
model  matching  results  are  extended  to  yield  a  solution  to  the  problem  of  tracking  by 
static  state  feedback. 


1  Introduction 

In  recent  years  there  has  been  a  great  deal  of  research  effort  in  the  adaptive  control  of 
nonlinear  systems.  This  research  has  been  primarily  focused  on  SISO  systems  for  which 
there  exist,  broadly  speaking,  three  types  of  approaches:  those  relying  on  the  existence  of 
certain  matching  or  structural  conditions  for  the  location  of  the  unknown  partuneters  (see 
for  example  [KKM89|,  [TKMK89]  and  [KKM91]),  the  second  relying  on  certain  assumptions 
on  the  type  of  the  nonlinearities  in  the  plant  (see  for  example,  [SI89],  [NA88],  [KTKS91]) 
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Abttract 

Stability  properties  of  zero  dynamics  are  among  the  crucial  input.output  proper* 
ties  of  both  linear  and  nonlinear  systems.  Unstable,  or  “non'minimum  phase”,  zero 
dynamics  are  a  major  obstacle  to  input«output  linearization  and  high  gain  designs.  An 
analysis  of  the  effects  of  regular  perturbations  in  system  equations  on  zero  dynamics 
shows  that,  whenever  a  perturbation  decreases  the  system's  relative  degree,  it  man¬ 
ifests  itself  as  a  singular  perturbation  of  zero  dynamics.  Conditions  are  given  under 
which  the  zero  dynamics  evolve  in  two  timescales  characteristic  of  a  standard  singular 
perturbation  form  that  allows  a  separate  analysis  of  slow  and  fast  parts  of  the  zero 
dynamics.  The  slow  part  is  shown  to  be  identical  to  the  zero  dynamics  of  the  unper¬ 
turbed  system,  while  the  fast  part,  represented  by  the  so  called  houndary  layer  system, 
describes  the  effects  of  perturbations. 
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Abstract 

An  indirect  adaptive  control  law  based  on  certainty  equivalence 
is  designed  for  a  model  of  the  induction  motor  with  the  as'^umption 
that  the  magnetic  subsystem  is  linear.  This  nonlinear  rcr.:rol  law  ren¬ 
ders  the  induction  motor  system  input-output  linear  and  also  achieves 
input-output  decoupling.  In  addition,  we  find  for  the  specific  case  of 
the  induction  motor  we  are  able  to  prove  parameter  convergence  and 
asymptotic  tracking  of  a  reference  trajectory  using  the  indirect  adap¬ 
tive  controUer.  This  result  differs  from  the  generic  case  where  we  can 
only  show  asymptotic  tracking.  The  indirect  adaptive  control  method¬ 
ology  also  does  not  suffer  from  the  drawback  of  overparameterization 
as  in  the  direct  adaptive  control  technique.  Simulations  are  also  given 
comparing  nonadaptive,  direct  adaptive,  and  indirect  adaptive  non¬ 
linear  controllers. 
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Abstract 

A  toolbox  for  nonlinear  control  system  desip  is  presented.  This 
package  contains  modules  to  approximate  systems  to  polynomials  sys¬ 
tems  of  arbitrary  order  and  then  render  them  input-output  linear  or 
input-state  linear  with  arbitrary  order  error  terms.  We  also  discuss 
possibilities  for  real-time  control. 
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In  this  paper  we  present  matching  conditions  for  output  regulation  for  four 
major  classes  of  perturbed  nonlinear  systems  controlled  via  the  method  of 
exact  linearization  utilizing  the  sliding  mode  control  methodology.  The  sys¬ 
tems  considered  are  single  input  single  output  (SISO)  systems  with  perturbed 
zero  dynamics,  multiple  input  multiple  output  (MIMO)  systems  with  well  de¬ 
fined  vector  relative  degree,  left  invertible  MIMO  systems  decoupled  using 
the  zero  dynamics  algorithm,  right  invertible  MIMO  systems  decoupled  using 
the  dynamic  extension  inethod. 
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